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I. INTRODUCTION 



This thesis supplements and updates ongoing research at the Naval Postgraduate School. 
It documents and describes work completed to improve existing software algorithms used to 
process holographic images of rocket motor combustion products. 

A. BACKGROUND 

Holographic techniques are used to produce an image of the rocket motor combustion 
products so that a statistical analysis of the particles in the gases may be performed. The 
particles are from additives to the solid propellant. They are added for the purpose of improving 
certain performance characteristics of the rocket motor. The particle size distribution has a 
major effect on the performance of the motor, and it is this information that is sought in the 
research. A complete description of the hologram recording, reconstruction, and imaging 
process can be found in References 1 , 2, and 3. 

The main objective of this thesis was to improve the existing microcomputer software 
programs used to perform automatic data retrieval from the holograms. The area of research 
was directed mainly at reducing the executable code size and improving the time of execution 
of the programs. Two phases of improvement were attempted. The first was the application 
of optimizing compiler techniques to the existing FORTRAN programs. The second was a 
subsequent complete rewrite of the programs using the C programming language and an 
optimizing compiler. 

B. HARDWARE DESCRIPTION 

The hardware system consists of a standard IBM PC/AT computer with 40MB hard drive 
configured with an INTEL INBOARD 386/AT microprocessor with a 80387 math coprocessor 
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module, an AST ADVANTAGE memory expansion board with 1.5 megabytes of memory, and 
an Imaging Technology Inc. PCVISION Frame Grabber image processing board installed in an 
eight bit data bus slot of the IBM PC/AT. A detailed description of the PCVISION Frame 
Grabber hardware can be found in Reference 4. The external support items include a video 
monitor and VCR recorder for input, viewing, and recording of external images; an external 5.25 
(20Mb) cartridge Bernouli disk drive used to store image files; and a Tektronix thermal copier 
used to get hardcopy of the displayed video monitor screen. 

C. SYSTEM SOFTWARE 

The software associated with the system includes the following: 

- Imaging Technology’s IMAGEACTION software package, a menu driven software package 
that provides numerous routines for real time processing of images [Ref. 5], 

- Imaging Technology’s ITEX/PC software library which contains the routines used in both the 

FORTRAN and C locally produced programs [Ref. 6]. 

- IBM Disk Operating System Version 3.3 (DOS). 

- Microsoft FORTRAN Compiler Version 3.3. 

- Microsoft FORTRAN Optimizing Compiler Version 4.1 (MSFORT) [Refs. 7, 8, and 9]. 

- Microsoft C Optimizing Compiler Version 5.0 (MSC) [Refs. 10, 11, and 12]. 

- Locally produced automatic data retrieval and analysis programs used to support the research 

[Refs. 1, 2, and 3]. 

D. THESIS ORGANIZATION 

This thesis consists of four chapters. Chapter Two discusses the improvements made in 
the software systems and also gives a detailed discussion of the new programs developed 
using the C language. Chapter Three discusses the software and hardware performance 
analysis and shows the results of improvements made. A summary and concluding remarks 
are presented in Chapter Four. Appendices A through J contain the source code for the C 
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Language programs developed. These programs are heavily documented to explain their 
operation. Appendix K contains a description and files used to simplify and automate FORTRAN 
and C compiler operations. 
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II. SOFTWARE IMPROVEMENTS 



Two phases of code optimization were conducted. The FORTRAN optimization was achieved 
using the MSFORT compiler version 4.1 [Ref. 7] and the C programming optimization was 
accomplished using MSC compiler versiort 5.0 [Ref. 10]. 

A. FORTRAN OPTIMIZATION 

Prior to this thesis, the Microsoft FORTRAN Compiler Version 3.3 was used exclusively to 
compile locally produced programs. Previous programming done by Redman [Ref. 1], Edwards 
[Ref. 2], and Orguc [Ref. 3] used this compiler to compile all FORTRAN programs. One of the 
results of this thesis was the elimination of the dependence on this non-optimizing compiler. 

The MSFORT optimizing compiler was available for use during the work done by Orguc but 
could not be used due to incompatibilities with the compiler and the ITEX/PC software library. 
The incompatibility was determined to be due to Microsoft’s updating several aspects of the 
compiler to conform with the FORTRAN 77 ANSI standard. One of the improvements to the new 
compiler involved stringent enforcement of type casting between different modules in a program 
(i.e., main and subroutine modules). This type checking requirement resulted in several 
problems in the original programs due to type casts of variable declarations (i.e., CHARACTER 
*21) that were not consistent across main program and subroutine modules (i.e., main program 
declaration of CHARACTER *21 and the subroutine declaration of CHARACTER *22). This 
resulted in errors when compiling the programs but was not very evident when cross referencing 
the error code produced by the compiler. These type casts did not produce compile errors 
when compiling with the old compiler version. 

The second problem dealt with errors generated within the ITEX/PC library FORTRAN include 
file itexpc.inc [Ref. 6] which is required to interface the frame grabber routines. The compiler 
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error code statement generated was F2115: syntax error. This error was attached to all array 
type declarations within the file. The error was traced to the format of the array declarations 
within the include file and was corrected by reversing the items within the declaration. A sample 
declaration before correction is 

integer*2 Jarray(9) [REFERENCE, NEAR], 
with the correction being 

integers )array[REFERENCE,NEAR] (9). 

Once all of the syntax errors due to the above errors were corrected, the MSFORT compiler was 
usable with the ITEX/PC Software package. 

B. C PROGRAMMING OPTIMIZATION 

The C language optimization consisted of a complete rewrite of the locally developed 
FORTRAN routines using the MSC compiler. Several major achievements were realized using 
the C language. One significant improvement was the elimination of the need to define large 
data arrays for image pixel value storage. These data arrays used all the available memory 
assets (640 Kb DOS Limit) of the computer and resulted in the programs having to be written 
to use several loops or repeats of code to accomplish analysis of an entire image (512x480 
pixels). In some cases the analysis of the entire image was not completed due to this 
restriction. The overall effects were mainly large-size executable programs that ran slower due 
to the extra overhead to compensate for these effects. 

The C language allowed for the use of dynamic memory allocation during runtime and 
repeated calls to the functions of the ITEX/PC C library without a significant degradation in 
program speed or size. The use of a disk-based virtual array program also allowed for the easy 
use and access of large data arrays, when required, with minimum effect on the stringent 
memory restrictions. As a result the programs now support analysis of the entire usable image 
(512x480 pixels) with no special programming efforts to accomplish the tasks desired. 
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The following paragraphs describe the C programs generated for this thesis research. The 
performance and statistical analysis of these programs are presented in Chapter 3. 

1 . Program File thesis.h 

The file thesis.h (Appendix A) is the header file for use with all the programs. It contains 
the required library include files necessary for successful compilation of the programs. It also 
serves to define the C manifest constants that are used throughout the programs. The main 
purpose for this header is to maintain a sense of portability among the constants and variable 
definitions used within the programs. It makes for easy redefinition of the constants if necessary 
without having to change these values wherever they appear within the code. 

2. Program File genfunc.c 

The general support functions program file genfunc.c (Appendix B) documents and lists 
a collection of functions used throughout the program files to perform routine evolutions. Each 
function title and comment within the appendix defines its purpose. Also included with this file 
is an array definition named p[x][y]. This array was used to test the programs offline before 
they were used on the image processing system. Use of the array is easy to implement and 
is explained in the comment statements located within the array definition. 

3. Program File threshit.c 

The Image threshold function program threshit.c (Appendix C) is the C language version 
of the threshold subroutine program presented by Edwards [Ref. 2]. This program takes an 
operator input value that is used as a threshold limit for the image and sets all pixel values 
below this value to BLACK (0) and those above it to WHITE (255). This effectively changes 
all pixel values representing image features to one value, thus producing a binary image of the 
original. This routine is used primarily to produce an image suitable for input to the feature 
identification program described later. 
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4. Program File speckle.c 

Speckle noise index function, speckle.c (Appendix D), is a routine that uses the 
algorithm presented in Reference 13 and used by Edwards [Ref. 2] to calculate the value of 
speckle index. The value is a measure of the speckle noise present in an image. It is used to 
evaluate the effectiveness of a filter operation that has been performed to reduce the amount 
of speckle noise present in the image. This routine is used by all the filter routines to serve as 
a measure of their effectiveness. 

5. Program File virjarry.c 

The virtual array support functions, file virjarry.c (Appendix E), are a collection of the 
functions necessary to setup and access the virtual disk-based array scheme used in the filter 
programs discussed later. These functions were originally presented in Reference 15 and 
required only minor modification to support this research effort. 

The access macro definitions and C type casts necessary to use the virtual array routines 
are defined in the program header file thesis. h (Appendix A). It is these definitions and macros 
that make the routines easy to use. When properly set up, the virtual disk array can be used 
as if it were the same as any C language array declared internal to a program file. All disk 
evolutions and file handling are done in the background without operator intervention. For a 
detailed analysis, or before modification of any of the functions, it is recommended that 
Reference 15 be consulted. 

a. Function init_vjarray (filename, record jsize, 

filchar) 

This function creates a virtual array file named filename on the default disk drive. 
The size of each element will be set to the value of record jsize. The variable filchar is the 
character that is set into initial array values before formal assignment of the value is made. 
The number of elements in the array is set initially to zero and the file closed. This routine 
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provides the necessary file and initial settings to be used by the other functions called to use 
the virtual array. For example, a statement such as 

l n H_ v _«rr«y C'plx. ya r,*l z eof (lr»t) ,N U LL) ; 

would create a file named pix.var that would have integer elements and use the defined NULL 
value as the fill character. 

b. Function open _v_array(filename, buffer jsize) 

This function prepares an existing virtual array for use by opening the existing file of 
filename. It dynamically allocates memory space for a buffer to hold the elements of the array 
for immediate use. The size of this buffer is the value bufferjsize. A statement such as 

array = open_v_arrayfpix.var“,100); 

would open the file pix.var for use and create a buffer to hold 100 elements. The variable 
array must have been defined as a pointer variable. 

c. Function close_v_array(array) 

The function close_v_array writes any remaining elements from the buffer to disk, 
closes the virtual array file, and releases the allocated memory for the buffer. This routine 
should be the last one called in the program before exit. 

d. Function access _y_rec(array, index) 

This routine performs the low-level access to the file and makes sure the array 
element referenced by index is available in the memory buffer. If the element is not in the 
buffer, it reads the element from the disk. It also will extend the size of the virtual array 
automatically if the indexed value is not in the array. This allows for the elements of the array 
to be assigned only when needed and the array file only to grow as large as necessary to 
hold declared array values. 

6. Program File featjd.c 

Program featjd.c (Appendix F) is a complete stand-alone program that is used to scan 
the image seeking connected features and to assign a unique identification number to each 
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identified feature within an image. The maximum number of features that can be counted is 
limited only by the maximum value that the integer variable can assume (32,767). This number 
is then used to give a count of the number of features present in the image. The program must 
be supplied an image that has been thresholded since all decisions within the program are 
based on the value of the pixels being either WHITE or BLACK. 

This program is a complete rewrite of the FORTRAN program presented by Orguc [Ref. 
3). It uses essentially the same program sequencing of events and ideas as the original, but 
the use of the GO TO statement that was used numerous times in the FORTRAN version to 
control program flow has been entirely eliminated from this version resulting in smoother 
program flow. 

7. Program File sizeit.c 

Image feature sizing is accomplished using program sizeit.c (Appendix G). This routine 
is a rewrite of the program presented by Orguc [Ref. 3]. It is a stand-alone program that takes 
the feature labeled output from the identification program ( featjd.c ) and produces an output 
file listing each feature’s area, x-chord, and y-chord. The size.dat output file is in a form 
suitable for input into a statistical analysis program such as STATGRAPHICS [Ref. 14]. This file 
must be saved under another filename or processed into STATGRAPHICS before the sizing 
program is executed again because any existing file is overwritten during program execution. 
The main new feature of the C language version over the FORTFIAN is that dynamic memory 
allocation was used to allocate data storage for each feature’s area, x-chord, and y-chord. This 
removed a restriction on the maximum number of features that could be processed. This 
program also uses the SCALE_FACTOR constant to convert the image pixel size data to values 
indicating the actual physical size of the particles. This constant value is based on the 
equipment set up during hologram processing. It is determined by making an image of a 
known object, such as a threaded screw, at the same magnification. This test image is then 
used to determine a calibration value for the conversion from pixels to actual feature size. 
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8. Program File Istat.c 

The local statistical filter algorithm Istat.c (Appendix H) is the first of three filters used in 
the analysis to reduce speckle noise in the images. This version is a C language version of the 
FORTRAN program presented by Edwards [Ref. 2]. It uses the disk virtual array functions [Ref. 
15] to provide storage for the filtered pixel array during program calculation. The algorithm is 
described in Reference 1 6 and basically uses a 5x5 array of pixel values around a central pixel 
and applies the statistical average obtained from this local array to the central pixel as a new 
pixel value. Once the calculation for all the pixels in the image is complete, the new pixel values 
are written to the image screen. 

9. Program File sigma.c 

The program sigma.c (Appendix I) is the C version of the 2sigma filter algorithm of 
Reference 16 and presented by Edwards [Ref. 2]. As with the local statistical filter, this program 
uses the virtual array functions to store the filtered pixel array during calculation in a disk based 
virtual array. The algorithm basically takes a 5x5 array average of pixel values and removes all 
values outside of the 2sigma range. This calculated value is then applied to the 5x5 array 
central pixel and stored for this central position. Once all pixel positions have been calculated, 
the new pixel values are written to the image screen. 

10. Program File geofil.c 

The program file geofil.c (Appendix J) is the C language adaption of the filter program 
introduced by Edwards [Ref 2]. This filter uses a geometric hulling algorithm [Ref. 13] to filter 
the image by comparing pixel values in the horizontal, vertical, and two diagonals around a 
central pixel. It compares the original image and its complement to determine the proper level 
to set the central pixel based on the neighboring pixel values. The overall result is the reduction 
of the speckle noise components at a faster rate than the actual feature pixel data. The 
complement image array is stored using the virtual disk array routine and thus fully supports 
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using the entire usable image (512x480). The original program only filtered the first 120 rows 
of the image (512x120) and provided no means for evaluating any other part of the image. 

C. SOFTWARE DEVELOPMENT TOOLS 

Several development tools were available for use during program development. Both the 
MSFORT and MSC compiler packages come with the Micro Soft CodeView Debugger program 
[Ref. 9 and 12]. This program is an online debugger that allows for program debugging 
provided special compiler directives are used when compiling the programs. 

Since use of CodeView was not easy to accomplish without the extra compiling, a special 
array p[x][y] (Appendix B) was developed to simulate the image screen. This array was then 
substituted where necessary to simulate use of the library routines that access the image pixels. 
With the substitution installed, the programs could be compiled and then run to check the 
results of program execution. A side effect of this process was the ability to also compile the 
programs using Borland’s Turbo C 2.0 [Ref. 17] compiler which includes an easy to use online 
source code debugger within an integrated environment. This process was easier to use than 
the MSC Codeview program. The TURBO C package is not part of the existing imaging system 
but should be considered for future upgrades to the system. The package was recently 
upgraded to version 2.0 and in some respects is better than the MSC package. 

Appendix K is a listing and explanation of the DOS Batch files used to compile and link the 
programs with each compiler. The compiler and linker command directives are those 
recommended in the ITEX/PC Programmer’s guide [Ref. 6]. An explanation of each directive 
is contained in References 7 and 10 respectively. 
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III. SOFTWARE PERFORMANCE 



This chapter deals specifically with the software performance (both FORTRAN and C) of 
the programs developed and improved during this research. The improvements in the 
FORTFtAN versions were the reduction of executable program size and the ability to run in the 
latest DOS environment (DOS Version 4.0). The original programs required the absence of 
any specially installed device drivers in the DOS environment and would not allow any TSR 
programs to be loaded during normal processing on the IBM PC/AT. The use of the optimizing 
compiler allowed all previous programs to execute while coexisting with what is considered a 
normal set of DOS installed drivers. The newid and newsize FORTRAN programs, however, 
still require all of the DOS memory limit for execution and must have any TSR programs 
removed from memory, but no change in the DOS environment is required before execution. 

The improvements using the C language were twofold. The first was a decrease in execution 
time, and the other was the ability to easily process the entire usable system image size 
(512x480 pixels) for the first time. The latter seems to be the most significant since future 
plans to upgrade to the latest ITEX/PC imaging hardware and software will allow an even larger 
image to be processed. 

No specific program testing or verification of the original FORTRAN programs was 
accomplished. This was deemed unnecessary since no actual coding changes were made in 
the programs. Only the compiler directives were changed to meet the new format introduced 
with this version of the compiler. 

The testing and verification of the C programs consisted of taking images that were 
previously processed, recreating a test environment to match that used for the original 
programs, and then comparing the results against the reported findings of Reference 2 and 
Reference 3. 
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A. PROGRAM SIZE INFORMATION 



The FORTRAN and C executable program size information is presented in Table 1 under the 
MSFORT and MSC headings. Each entry is for the EXE (executable file) listing as obtained 
from the DOS DIR command. The size given is in Kilobytes (KB) for each file as it is stored on 
disk. The missing table entry for the MSC version of SPECKLE is because this routine is being 
used in each filter and is included in each C program filter. 

The dramatic decrease in size from the original FORTRAN versions to those of the MSFORT 
compiler are believed to be mostly attributable to the use of the EXEPACK utility of the new 
compiler [Ref. 7]. The sizes for the MSFORT and MSC versions are relatively consistent with 
the exception of the filter programs. The added speckle index algorithm and the virtual array 
code in each filter accounts for the C filters being larger than their FORTRAN counterparts. 



TABLE 1 

EXECUTABLE PROGRAM SIZE 



Filename 


Original 


MSFORT 


MSC 


THRESH . EXE 


131,882 


44,795 


31,467 


SPECKLE . EXE 


122,180 


33,569 


NONE 


NEWSIZE.EXE 


556,038 


52,575 


46,281 


NEWID.EXE 


569,830 


46,655 


32,509 


SIGMA. EXE 


379,704 


48,293 


49,483 


STAT.EXE 


544,626 


44,143 


49,475 


GEOFIL.EXE 


364,382 


41,375 


49,969 
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B. C PROGRAM VERIFICATION AND EXECUTION TIME COMPARISON 



The testing and verification of the new C programs consisted of using existing images that 
were previously analyzed and reported on. The two images chosen are shown in Figure 1 and 
Figure 2. The image lOxwgrid.img (Figure 1) was chosen from the image library to test the 
operation of the featjd.c and sizejt.c programs. This image shows 23 features that are taken 
from the standard calibration image used by Redman [Ref. 1]. It displays the 23 different 
feature sizes that are present in the calibration image. Figure 2 is the j17res4.img library image 
and is an image taken of the AFRT with speckle noise present. This image was used to test 
each filter program for proper operation. In each test the results obtained were compared 
against data previously obtained for this project. 

The lOxwgrid.img image was thresholded at a value of 130 and then processed with both 
the FORTRAN and C versions of the sizing and identification programs. The results of these 
tests appear to indicate proper operation of the new programs. Both featjd.c and sizejtc 
produced the same outputs as programs newid.exe and newsize.exe [Ref. 3]. Both counted 
and sized 23 features. Tables 2 and 3 show the data output for the two programs. The results 
were normalized to remove the physical size scaling and calculations for conversion to microns. 
Normalization was required since there does exist a discrepancy between how the programs 
calculate the final physical size data. The size data obtained from sizejt.c is about a factor 
of four smaller than that obtained with the program newsize. The SCALE_FACTOR entry in the 
thesis.h header file or the magnification calculation in newsize may be incorrect. These values 
should be adjusted or corrected before a quantitative value for the feature size is used in any 
critical calculations. 
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Figure! Library Image 1 0xwgrid.img 




Figure 2 Library Image j17res4.img 
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TABLE 2 

PROGRAM NEWSIZE OUTPUT DATA. 



Num 

I 


Area 

1 


X-chord 

l 


Y-chord 


1 


101.002 


13.000 


4.000 


2 


10.000 


5.000 


2.000 


3 


3.000 


2.000 


1.000 


4 


6.000 


3.000 


2.000 


5 


15.000 


4.000 


4.000 


6 


20.000 


5.000 


6.000 


7 


70.002 


9.000 


10.000 


8 


106.003 


11.000 


12.000 


9 


138.003 


12.000 


14.000 


10 


177.004 


14.000 


16.000 


11 


199.005 


15.000 


17.000 


12 


262.006 


17.000 


20.000 


13 


310.007 


19.000 


21.000 


14 


377.009 


21.000 


24.000 


15 


429.010 


22.000 


25.000 


16 


478.011 


23.000 


27.000 


17 


562.013 


25.000 


29.000 


18 


626.015 


26.000 


30.000 


19 


733.017 


29.000 


33.000 


20 


861.020 


31.000 


35.000 


21 


948.023 


31.000 


39.001 


22 


1200.029 


36.000 


43.001 


23 


1468.035 


40.001 


46.001 
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TABLE 3 

PROGRAM SIZE IT OUTPUT DATA. 



NUM 

l 


AREA 

i 


X-chord 

i 


Y-chord 


1 


1468.000 


40.000 


46.000 


2 


1200.000 


36.000 


43.000 


3 


948.000 


31.000 


39.000 


4 


861.000 


31.000 


35.000 


5 


733.000 


29.000 


33.000 


6 


626.000 


26.000 


30.000 


7 


562.000 


25.000 


29.000 


8 


478.000 


23.000 


27.000 


9 


429.000 


22.000 


25.000 


10 


377.000 


21.000 


24.000 


11 


310.000 


19.000 


21.000 


12 


262.000 


17.000 


20.000 


13 


199.000 


15.000 


17.000 


14 


177.000 


14.000 


16.000 


15 


138.000 


12.000 


14.000 


16 


106.000 


11.000 


12.000 


17 


70.000 


9.000 


10.000 


18 


20.000 


5.000 


6.000 


19 


15.000 


4.000 


4.000 


20 


6.000 


3.000 


2.000 


21 


2.000 


2.000 


1.000 


22 


10.000 


5.000 


2.000 


23 


1.000 


1.000 


1.000 
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The j17res4.img image was used to test each filter program. The results of the testing are 
presented in Table 4, Table 5, and Figure 3. Each program was run using a DOS text input 
file that used DOS redirection to supply the keyboard inputs to the program without operator 
input. An external software timer measured the times for loading and executing the programs. 
The timing data in Table 4 were all gathered in this manner. Table 5 is a tabulation of Speckle 
Index values obtained during each filter operation and contains the data plotted in Figure 3. 
These data clearly indicate agreement with that achieved previously by Edwards [Ref. 2] and 
Orguc [Ref. 3]. 

Figures 4 and 5 show one iteration of the 2sigma filter applied to a full screen image 
(512x480). The image is once again the j17res4.img file that has a pixel value histogram 
generated by the IMAGEACTION package imbedded into the picture. Comparison of the two 
histograms visually shows that the filter reduces the speckle content of the image as expected. 
The presence of the valley and hump to the left of the histogram is indication that the feature 
grey scale data is emerging above the speckle noise. This result agrees with that achieved by 
Edwards [Ref. 2]. 



TABLE 4 

PROGRAM EXECUTION TIME COMPARISON 



Program 


MSFORT 


MSC 


Feat ID 


1 m 11 s 


1 m Os 


Feat Sizing 


1 m 15 s 


1 m 35 s 


2Sigma Filter 


7 m 20 s* 


6 m 3 3 s 


Local STAT Filter 


39 m 45 s* 


13 m Os 


Geometric Filter 


18 m 50 s* 


6 m 17 s 



* [Ref. 2] 
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TABLE 5 

SPECKLE INDEX REDUCTION RESULTS 



NUM 


2 SIGMA 


L STAT. 


GEOMETRIC 


0 


0.304005 


0.304005 


0.304005 


1 


0.153767 


0.194121 


0.224875 


2 


0.098967 


0.151951 


0.178222 


3 


0.074415 


0.129769 


0.146152 


4 


0.060615 


0.115938 


0.122347 


5 


0.051691 


0.106372 


0.103983 


6 


0.045354 


0.099367 


0.089425 
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Figure 4 Unfiltered Image jl 7res4.img 
(SI = 0.304005) 




Figure 5 Filtered Image j17res4.img 
(SI = 0.153767) 
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IV. CONCLUSIONS 



In summary, this thesis has achieved a significant increase in the ability to perform 
hologram image analysis on the present imaging system. The goals of smaller compact 
programs coupled with faster execution times were both achieved. The capability to perform 
analysis of full 512x480 images in all programs is now available. 

Another important aspect is that of portability to the recently acquired upgraded system that 
includes a COMPAQ 386/20 computer and a new IMAGING Technology Processor and software 
that only uses the C Language. This new system will be the subject of future research and 
these programs should port with minimum effort over to the new system. The only requirement 
should be verification of the required function calls for the library routines used with those of 
the new library and then the changes made to the header file thesis.h. No major program 
code changes should be required. 

It is recommended that further verification testing be conducted with these programs using 
other images. The testing conducted for this thesis only showed that the programs operate 
properly and produce results similar to those achieved previously. 
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APPENDIX A 



PROGRAM HEADER FILE: thesls.h 



/* The purpose of this file is to list all necessary include files, 
manifest constants, and MACRO definitions required by all the 
programs within this package. The programs were written to be as 
portable as possible and, as such, changes to constants in 
this file should be all that is required to change a parameter 
definition throughout the program files. *1 

f* Include files for use with ITEXPC programs */ 



#include <stdio.h> 

#include < stdlib.h > 

#include <math.h> 

#include <malloc.h> 

#include 'Hexpc.h' 

*/ 

/* Hex board base memory start address */ 
/* Hex board base register start address */ 

/* Hex board color option installed */ 



/* InHial ITEX/PC Board Jumper Settings 

#define MEMBASE OxDOOO 
#define REGBASE OxFFOO 
#define IFU\G PSEUDO_COLOR 



I* ITEXPC inHial AOI (area of interest) settings */ 



#define IXS 0 

#define IYS 0 

#define NROW 480 
#define NCOL 512 



f* HexPC LUT (Look Up Table) Variables */ 

f* color board variables */ 

#define RED 0 

#define GREEN 1 
#define BLUE 2 

#define INPUT 3 



f* Threshold lim'rts */ 



#define 


LOWEST 


0 


/* Equates to Black for lowcut value */ 


#define 


HIGHEST 


255 


1* Equates to WhHe for highcut value */ 


#define 


BLACK 


0 




#define 


WHITE 


255 




/* FiHer array sizes 


and limHs */ 




#define 


SNUM 


25 




#define 


NUM 


9 




#define 


HIGH 


254 




#define 


LOW 


0 





/* Sizing magnification factor used in sizing program */ 

#define SCALE FACTOR 3.7353 /* Based on image calibration V 
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/* Virtual array Header File Definitions */ 
#define header 7 



/* Virtual Array Control Block typedef */ 
typedef struct { 



FILE Mile; 
long size; 
int elsize; 
char ^buffer; 
int buf_elsize; 
int buf_size; 
char *blank_rec; 



} 



/* pointer to file control block */ 

/* number or array elements in file */ 

/* number or bytes in each element V 
I* pointer to array buffer */ 

/* size of element in buffer including index V 
/* number of elements in buffer */ 

/* pointer to initialization record */ 

I* used for extending file V 



VACB ; /* Virtual Array control block type name V 



t* Virtual Array Access Prototypes *1 



int init_v_array(char Milename.int rec_size,char filchar); 
VACB *open_v_array(char Milename.int buffer_size); 
void close_v_array(VACB *v_array); 
void *access_v_rec(VACB *v_array,long index); 

/* Virtual Array Access Macros */ 

#define VREC(y) ((Kerns *)access_v_rec(item_array,y)) 
#define gdesc(y) VREC(y)->v_gdesc 
#define G(x,y) VREC(y)->v_gdesc[x] 

/* Virtual Array element structure typedef V 

typedef struct { 

unsigned char v_gdesc[NCOL]; 

} 

items; 



/* 



*/ 



23 



APPENDIX B 



PROGRAM FILE: genfunc.c 



/”** 

TITLE : Image Processing General Support Functions 

CALLED BY: various routines 

FILENAME : genfunc.c 



UST 

MODIFIED : 12/8/1988 



PURPOSE : 

These general functions provide for routine evolutions 
that occur a number of times in tne Image processing 
routines developed for analysis of the speckle reduction 
algorithms. Some require support of other functions and 
are not totally independent. This file serves as 
documentation only and has not been nor can it be compiled 
as a stand-alone program file. 
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/* initial ITEXPC Board setup *1 

startitO 

{ 

sethdw(REGBASE,MEMBASE.IFLAG) ; 
aclear(IXS,IYS,NCOL,NROW,1 50) ; 
initializef); 

system (’els'); /* Clears monitor screen *1 

} 

/* Get integer keyboard input */ 

geti(name.iptr) 
char *name; 
int *iptr; 

{ 

printff %s '.name); 
scanff %d',iptr); 

} 

/* Get floating point keyboard input */ 

getf(rname.riptr) 
char *rname; 
int *riptr; 

{ 

printf(“%s ",rname); 
scanff %f,riptr); 

} 

f* Pause routine used to prompt for operator intervention*/ 
wait() 

{ 

printff\n\n Press Return to continue 1 ); 
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fflush(stdin); 

getcharO; 

} 

f* returns minimum integer value from input array passed */ 

getmin (data, num) 

int *data; 
int num; 

{ 

int t, min; 

for (min= data[0], t=1 ; t<num ; t++) 
if( data[t] < min ) min = data[t]; 

return min; 

> 

/* returns maximum integer value from input array passed */ 

getmax(data,num) 

int *data; 
int num; 

{ 

int t, max; 

for (max= data[0], t=1 ; tcnum ; t++) 
if( data[t] > max ) max = data[t]; 

return max; 

} 

/* DOS System Call to clear Display Screen */ 

clsO 

{ 

systemf'cls"); 

> 

I* */ 

/* end of function programs *f 



I* TESTER */ 

/* The following defined array was and can be used to test routines 
before they are applied to the imaging system (i.e., totally offline 
with any computer). This particular array was used in the debugger 
program to verify proper operation of the algorithms before they were 
tested using the full image. To use the tester, include it in the 
program and change the calls to RPIXEL and WPIXEL as needed. 



Change rpixel(x,y) to p[y][x] and use the following data array 

NOTE: The order of subscript changes ’[y][x]’ is due to the way C calls 
an array. 

To simulate different screen attributes, change array values as 
needed to simulate the image screen desired. */ 
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static int p[15] 



{ 1. 2. 3. 4 

( 1. 2. 3, 4 

C 1. 2. 3, 4 

{ 1. 2. 3, 4 

c 1. 2, 3, 4 
€ 1. 2, 3, 4 

< 0 , 0 , 0 , 0 

< 255, 4, 4,255 

< 255, 4, 4, 4 

{ 255,255,255, 4 

< 255,255,255, 4 

< 255,255, 4, 4 

{ 255,255,255,255 
C 255,255,255,255 
{ 255,255,255,255 
>; 



[ 20 ] = { 

, 5, 0, 

, 5, 0, 

, 5, 0, 

, 5, 0, 

, 5, 0, 

, 5, 0, 

, 0 , 0 , 
,255,255, 
, 4, 4, 

, 4, 4, 

, 4,255, 

, 4, 4, 

,255,255, 
,255,255, 
,255,255, 



255,255 

255,255 

3,255 

3,255 

3,255 

255,255 

255,255 

255,255 

255,255 

255,255 

255,255 

255,255 

255,255 

255,255 

255,255 



,255,255, 
,255,255, 
,255,255, 
,255,255, 
,255,255, 
,255,255, 
,255,255, 
,255, 5, 

5, 5, 
,255, 5, 

,255,255, 
,255,255, 
,255,255, 
,255,255, 
,255,255, 



255,255 
255,255 
255,255 
255,255 
255,255 
255, 5 

5 , 

5 , 

5 , 

5 , 

5 , 
255, 
255,255 
255,255 
255,255 



, 255, 255 
, 255, 255 
,255,255 
,255,255 
,255,255 
, 5,255 

- 5, 5 

- 5, 

, 5, 

- 5, 

- 5, 

, 5,255 

,255,255 
,255,255 
,255,255 



,255,255 
,255,255 
,255,255 
,255,255 
,255,255 
,255,255 
,255,255 
5,255 
5, 5 

5,255 
,255,255 
,255,255 
,255,255 
,255,255 
,255,255 



,255,255, 
,255,255, 
,255,255, 
,255,255, 
,255,255, 
, 6 , 6 , 
, 6 , 6 , 
, 6 , 6 , 
, 6 , 6 , 
,255,255, 
,255,255, 
,255,255, 
,255,255, 
,255,255, 
,255,255, 



255,255}, 
255,255}, 
255,255}, 
255,255}, 
255,255}, 
6 , 6 }, 
6 , 6 }, 
6 , 6 }, 
6 , 6 }, 
255,255}, 
255,255}, 
255,255}, 
255,255}, 
7, 7}, 

255, 7} 



/* 



END TESTER 



*/ 
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APPENDIX C 



PROGRAM FILE: threshtt.c 



I* *************************************************** *********************** 

TITLE : Image THRESHOLD Routine 
FILENAME : threshit.c 



LAST 

MODIFIED : 12/9/88 



PURPOSE : 

This program thresholds the image on screen by taking 
an operator input value and forcing all image pixel 
values above the threshold value to the value of BLACK 
and all those below the value to the WHITE value. 

( Background == WHITE ; Feature == BLACK) 



********************************************************************** 
#include Ihesis.h* 

LS20 filename; /* LS20 & LS200 are required for ItexPC */ 

LS200 comline; /* readimO and saveimO library routines */ 

mainO 

{ 

int flag; 

printff\nThis program will threshold the input image desired. 1 ); 
printff'VAnWeady to Load Image?.. .Yes (y) / No (n) '); 
flag=getchO; 

if (flag == T|| flag == ’y*) 

{ 

startrtO ; 
readitO; 

> 

if (flag = = 'N' 1 1 flag = = ’n 1 ) 

{ 

initialize)); 

exit(0); 

} 

printf("\n\n\tTHRESHOLD the image?...Ves (y) / No (n) ■); 
flag « getchO; 

rf(flag = = Y|| flag = = 'y 1 ) threshrtO; 

} 

/* */ 

threshitO 

{ 

int ans, c, option; 
unsigned int val; 
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while(l) 

{ 

initialize!); 

geti("\n\n\tENTER NEW THRESHOLD LEVEL (0-255) :\&val); 

setlut(INPUT.O); 

threshold (val, HIGHEST); 

setlut(RED.O); 

threshold(val, HIGHEST); 

setlut(GREEN f 0); 

threshold(val, HIGHEST); 

setlut(BLUE,0); 

threshold (val, HIGHEST); 

setlut(INPUT,0); 

printff\n\nSATISFIED WITH THRESHOLD AT %d ?...Yes (y) / No (n)...\val); 
ans = getchO; 

if(ans == T 1 1 ans == V") break; 

} 

printff\n\t\t 0: SAVE THE THRESHOLDED IMAGE ?..."); 

printf (*\n\t\t 1: QUIT (Leave thresholded Image for further processing.)..."); 

printf (*\n\t\t 2: RESTORE System to original input image..."); 

geti(*\n\t Select option by NUMBER: \&option); 

switch(option) 

< 

case 0: 

maplut(IXS,fYS,NCOL,NROW); /* Map the image to frame memory */ 
save it 0 ; 
break; 
case 1 : 

maplut(IXS,IYS,NCOL,NROW); f* Map the image to frame memory */ 
break; 
case 2 : 

initializeO; I* initialize without a mapping operation */ 
break; 
default: 
break; 



} 

} 

/* V 

startitO /* initial fTEXPC Board setup */ 

{ 



sethdw(REGBASE,MEMBASE,IFLAG); 
aclearflXS.IYS.NCOL.NROW, 1 50) ; 
initializeO; 

system("cls"); /* Clears monitor screen */ 

> 

i* */ 

readhO 

{ 

int errval; 
int ch; 

while (1) 

{ 

printf("\n\n ENTER IMAGE FILENAME as (DEV:name.lMG)->"); 

fflush(stdin); 

gets(filename.s); 

filename. I = (char)strlen(filename.s); 

errval = readim(IXS,IYS,NCOL, NROW, filename, &comline); 
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iff errval = = 0 ){ 

printff\nFILE-> %s \n\nCOMMENT: %s ', filenames, comline.s); 
break; 

> 

if ( errval ! = 0 ) 

{ 

printffVAaError reading Filel I*); 
iff errval == -2 ) prlntff\n\tFile Not Foundl“); 
iff errval == -3) printf('\n\t\tBad File Formatl*); 
printffTry Again?? ...Yes (y) / No (n) *); 
ch = getchf); 

iff ch = = T 1 1 ch =- V’ ) continue; 
iff ch == 'N’ 1 1 ch == ’n’ ) exrtfO); 

} 

} 

} 

I* */ 

savehf) 

{ 

int comflag, errval; 
int ch; 

initializef); 

whilefl) 

{ 

printf(*\n\n ENTER FILENAME... ( DEV: name.lMG )->'); 

fflushfstdin); 

getsffilename.s); 

filename.l = (char)strlen(filename.s); 

getif\n Enter File Compression value ?...( 0/1 ) ',&comflag); 

printff \n\n Enter Image Comment or Return ..to continue.. (Max 200 CHAR):...'); 

fflushfstdin); 

getsfcomline.s); 

comline.l = (char) strlenfcomline.s); 

if (comflag == 1) printff\n\nStoring Image using Compression - please wart!!'); 
if (comflag == 0) printff\nStoring Image - plese wait!! “); 

errval = saveimf IXS,IYS,NCOL,NROW, comflag, filename, comline ); 

iff errval == 0) 

{ 

printff\nlmage Save completed SAT..'); 
break; 

} 

iff errval 1=0) 

{ 

printf(*\n\a Error saving File!!'); 

iff errval == -1 ) printff ^Insufficient Disk space' ); 

printff\n\tTry Again?? ...Yes (y) / No (n) “); 

ch = getchf); 

rff ch == V || ch == V’ ) continue; 
iff ch == ’N’ || ch == ’n’ ) break; 

} 

} 

} 

/* V 

geti(name,iptr) /* Read/Get integer keyboard input */ 
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char *name; 
int *iptr; 

{ 

printff %s .’Vname); 
scanf(" %d',iptr); 

> 

/* */ 
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APPENDIX D 



PROGRAM FILE: speckle.c 

jit ******************************************************* ******************* 

TITLE : Speckle Noise Reduction Function 
FILENAME : speckle.c 
CALLED BY; All Filter routines 



LAST 

MODIFIED : 11/27/1988 



PURPOSE : 

Provides for calculating the speckle index of an image. 

The resulting value is used as a measure to evaluate the 
success of filtering an image to reduce the speckle noise. 

This routine is called in all filter algorithms. 

-> OPERATOR must enter a Number of ROWS and COLUMNS 
THE PROGRAM uses to calculate the speckle index. 

ROWS — > Maximum = 480 Minimum = 1 
COLS — > Maximum = 512 Minimum = 1 

************************************************************************** *j 



speckleO 

{ 

extern int srow, scol; 
int m, n, x, y, nn; 
long smax, smin; 
long sdata[NUM]; 

unsigned long deviation, ssum; 
float smean, slocal, stotal, tot; 
float spklindex; 

if( srow >= NROW ) srow = NROW; 
if( scol >= NCOL ) scol = NCOL; 

tot = (long)srow * scol; 

I* Commence calculation */ 

printf("\n\n\tCalculating Speckle INDEX..."); 

for ( n = 1 ; n < srow-2 ; n++) 

{ 

for ( m = 1 ; m < scol-2 ; m++) 

{ 

smax = 0; 
smin = 260; 
ssum = 0; 
nn = 0; 
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for( y = ( n-1 ) ; y < (n+2) ; y++ ) 

{ 

for( x = ( m-1 ) ; x < ( m+2 ) ; x++) 

{ 

sdata[nn] = rpixel( x,y ); 

ssum + = (long)sdata[nn]; 

rf( smax < sdata[nn] ) smax = sdata[nn]; 

iff smin > sdata[nn] ) smin = sdata[nn]; 

nn+ + ; 

} 

} 

deviation = smax - smin; 
iff ssum == 0 ) smean = 1; 
smean = ssum / NUM ; 
slocal = (float) deviation / smean; 
stotal + = slocal; 

> 

} 

spklindex = stotal / tot; 

printff\n\n\tThe Calculated %d by %d speckle index is %f '.srow.scol, spklindex) 



/* END Speckle Function *1 
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APPENDIX E 



PROGRAM FILE: vlr_arry.c 

y* ********** **************************************************************** 

TfTLE : Virtual Array Support Functions 
FILENAME : vir_arry.c 



LAST 

MODIFIED : 11/27/1988 



PURPOSE : 

Provides for the setup of a disk drive virtual 
array that can be indexed as if it were in the calling 
program’s data storage area. The routines provide for disk 
file access when required to retrieve data elements from the 
array. This allows the size of a declared array to be limited 
only by the amount of DISK SPACE available. See thesis.h' 
header file for setup and MACRO routines that support these 
functions. Also, program Istat.c fully implements the array 
routines. 

******************************************************************** *******y 

#include thesis.h' 



/* Virtual Array Access Routines */ 

int init_v_array (filename, recsize.filchar) 
char ^filename, filchar; 
int rec size; 



{ 



long size; 

FILE M; 

f = fopen(filename p , wb , ); 
if (f != NULL) 

{ 

size = 0; 
fwnte(&$ize,4,1,f); 
fwrite(&rec_size,2,1 ,f); 
fwrhe(&filchar,1,1,f); 
fclose(f) ; 
return (1); 

> 

else 



/* write array size of zero *f 
/* and array element size */ 
/* and fill char */ 

/* to file header *f 



return (NULL); 

> 

VACB *open_v_array (filename, buffer_size) 
char Milename; 
int buffer_size; 

{ 

VACB *v_array; 
char *buf_ptr; 
int i; 

char filchar; 
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r allocate virtual array control block */ 

varray = (VACB *) malloc(sizeof(VACB)); 
if (v_array == NULL) return(NULL); 

/* open virtual array file */ 

v_array->file = fopenffilename/r+b*); 
if (v_array->file == NULL) 

{ 

free(v_array); 

return(NULL); 

}; 



I* get array size and element size for control block */ 

f read (&v_array- > size, 4 , 1 ,v_array- >f ile) ; 
f read (&v_array- > el size, 2, 1 ,v_array- > f i le) ; 
fread(&filchar,1 ,1 ,v_array->file); 
v_array->buf_elsize = v_array->elsize + 4; 

/* allocate buffer */ 



v_arr ay- > buffer = (char *) malloc(v_array->buf_elsize * (buffer_size + 1)); 
if (v_array-> buffer = = NULL) 

{ 

fclose(v_array->file) ; 
free(v_array); 
return (NULL); 

}: 

v_array->buf_size = buffer_size; 

/* set up blank_rec using the fill character in array header */ 

/* for initializing new array elements */ 

buf_ptr = v_array-> buffer + v_array->buf_elsize * v_array->buf_size; 

v_array->blank_rec = buf_ptr + 4; 

for (i = 0; i < v_array->buf_elsize; i++) 

*buf_ptr+ + = filchar; 

/* set record index negative for all buffer elements V 

buf_ptr = v_array-> buffer; 

for (i = 0; i < v_array->buf_size; i++) 

{ 

♦((long *)buf_ptr) = -1L; 
buf_ptr + = v_array->buf_elsize; 

}; 

return (v_array); 



void close_v_array(v_array) 
VACB *v_array; 

{ 

int i; 

char *buf_ptr; 

long recjndex, file_offset; 

buf_ptr = v_array-> buffer; 



/* flush buffer */ 
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for (i = 0; i < v_array->buf_size; i++) 

{ 

/* check each element index *f 

I* if element present; write it to disk */ 

recjndex = *((long *)buf_ptr); 
if (recjndex >= 0) 

{ 

file_offset = header + recjndex * v_array->el$ize; 
fseek (v_array->file,file_offset,0) ; 

(write (buf_ptr + 4, v_array->elsize,1, v_array->file); 

>; 

buf_ptr + = v_array->buf_elsize; 



void *access_vjec(v_array, index) 
VACB *v_array; 
long index; 

t 

char *buf_ptr; 
int bufjndex; 

long recjndex, tempjndex; 



/* calculate buffer address of referenced element */ 

bufjndex = index % v_array->buf_size; 

buf_ptr = v_array*> buffer + bufjndex * v_array->buf_elsize; 

recjndex = *(long *)buf_ptr; 

/* if element present, return its buffer address V 
if (recjndex == index) return (buf_ptr + 4); 

I* if element doesn’t exist, extend the file *1 

if (index >= v_array->size) 

{ 

f see k(v_array»> file, 0,2); 

for (tempjndex = v_array->size; tempjndex++ <= index; ) 
fwrhe(v_array- > blank rec, v_array->elsize, 1, v_array->file); 
v_array->size = index + 1; 
fseek (v_ar ray- > f i le,0,0) ; 
fwrite(&v_array->size, 4, 1, v_ar ray- > file); 

}; 

I* if buffer slot is occupied by another element, V 
r save it to disk */ 

if (recjndex >= 0) 

{ 

fseek (v_arr a y- > file, recjndex * v_array->elsize + header, 0); 
fwrite(buf_ptr + 4, v_array->elsize, 1, v_array->file); 

}; 

/* read referenced element into buffer slot V 

fseek (v_array->file, index * v_array->elsize + header, 0); 



}; 

free (v_ar ray- > buffer) ; 
fcl ose (v_ar ray- > f i I e) ; 
free (v_ar ray) ; 



/* de-al locate buffer V 
f* close array file V 
f* de-allocate VACB */ 



> 
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fread(buf_ptr + 4, v_array->elsize, 1, v_array->file); 
*((long *)buf_ptr) = index; 

/* return address of element in buffer */ 

return (buf_ptr + 4); 

} 

/* */ 
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APPENDIX F 



PROGRAM FILE: featjd.c 

I * ***************** ******************** ****** ******************************* 

TITLE : Image FEATURE Identification Algorithm 
FILENAME : featjd.c 



LAST 

MODIFIED : 12/9/1988 



PURPOSE : 

Labels and identifies each feature in an image on the ITEX 
system. Reads pixel-by-pixel and groups the objects by 
assigning a unique ID number to each feature or object so 
that they can be processed by other routines. 

The ID routine requires a thresholded image on screen or 
input of a previously saved thresholded image from disk or 
input from disk and then thresholding. This program module 
includes an optional call to threshold if desired. 

( Background == WHITE ; Feature == BLACK) 

*********************************************************************** ****y 

#include "thesis. h' 

LS20 filename; I* LS20 & LS200 are required for ItexPC *J 
LS200 comline; /* readimQ and saveimQ library routines */ 



mainO 

int flag; 

printf(Vi\aYou MUST USE A THRESHOLDED IMAGE for this Program!!"); 
printf(“\n\n\tReady to Load image from disk?. ..Yes (y) / No (n) ’); 
flag=getch(); 

if (flag == T|| flag == V) 

{ 

start itO; 
readitO; 

> 

if(flag == 'N' 1 1 flag = = ’n’ ) ex'rt(O); 



printff\n\n\tNeed to THRESHOLD the image?...Yes (y) / No (n) ■); 
flag = getchO; 

if (flag == T|| flag == Y) threshitO; 
feat_idO; 

printffVAnSave image to Disk File?..Yes (y) / No (n) ■); 
flag=getchO; 

if (flag — 'V'l | flag == V) saveitQ; 
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7 



featjdf) 

{ 

register int x,y; 

int a,x1 f y1 ,x2,y2,x3,y3,x4,y4; 

int fid = 0; 

int gid = 0; 

int temp_fid = 0; 

int temp gid = 0; 

int maxfl.nl ,n2,n1 a.nl b,n1 c.kk.kkl .nnl ; 

printff\n\n\tSTEP ONE in Progress..."); 

for( y = 1 ; y < NROW ; y++ ) 

{ 

for( x = 0 ; x < NCOL ; x++ ) 

{ 

iff rpixelfx.y) == WHITE) continue; 
iffx == 0) 

{ 

wpixelf x.y.fid ); 
continue; 

} 

if( ( rpixel( x-1,y ) != WHITE) ) 

{ 

wpixelf x, y,r pixel (x-1,y) ); 
continue; 

} 

iff rpixelfx.y-1) != WHITE ) 

{ 

wpixelf x.y.rpixelf x,y-1 ) ); 
continue; 

> 

wpixelf x.y.fid ); 
a = x + 1; 



while (1) 

{ 

iff rpixelf a.y ) == WHITE ) 

{ 

wpixelf x.y.fid ); 
stepf &fid,&gid ); 
break; 

} 

iff rpixelfa.y-1) != WHITE ) 

{ 

wpixelf x.y.rpixelf a,y-1 ) ); 
break; 

} 

a + = 1; 

} 

} I* end x - for V 
} I* end y - for */ 

fid += gid*HIGH; 

printf("\n\n\tSTART STEP TWO... "); 

printf(^n\nEnter the BEST estimate of the max feature length"); 
getifVAnlt’s Better to be too large than too small! !\&maxfl); 
printfOn\n\tRUNNING"); 

forf y = 1 ; y < NROW ; y++ ) 
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for( x = 1 ; x < NCOL ; x++ ) 

{ 

nl = rpixel( x,y ); 
n2 = rpixel( x,y-1 ); 

rf( ( nl 1 1 n2 == WHITE ) 1 1 (nl == n2 )) continue; 



x3 = x - maxfl; 
x4 = x + maxfl; 
y3 = y - maxfl; 
y4 = y + maxfl; 

if( x3 < LOW) x3 = LOW; 
rf( x4 > NCOL ) x4 = NCOL; 
if( y3 < LOW) y3 = LOW; 
rf( y4 > NROW ) y4 = NROW; 

for ( yl = y3 ; yl < y4 ; yl + + ) 

{ 

for ( xl ■ x3 ; xl < x4 ; x1 + +) 

{ 

nla = rpixel( xl.yl ); 

if ( nla == WHITE) continue; 

if ( nla = = nl ) wpixel( x1,y1,n2 ); 

> 

fid -= 1 ; 

> 

> 

> 

printf("\n\n\tSTEP THREE in Progress'); 

for ( y = 1 ; y < NROW ; y+ + ) 

{ 

for( x = 1 ; x < NCOL ; x++ ) 

{ 

kk = rpixel( x,y ); 

if( kk == WHITE ) continue; 

kkl = kk; 

nnl = tempjid + temp_gid * HIGH; 
if( kkl < nnl ) continue; 
if( (kkl - nnl) < -350 ) kkl = kk + 508; 
if( (kkl - nnl) >= -350 && (kkl - nnl) < -150 ) 
kkl = kk + 254; 
if( kkl != nnl ) 

{ 

x3 = x - maxfl; 
x4 = x + maxfl; 
y4 = y + maxfl; 

rf( x3 < 0) x3 = LOW; 

if( x4 > NCOL-1 ) x4 = NCOL -1; 

if( y4 > NROW-1 ) y4 = NROW -1; 

for( yl = y ; yl < y4 ; y1 + + ) 

{ 

for( xl = x3 ; xl < x4 ; x1 + +) 

{ 

nlc = rpixel( x1,y1 ); 

if( ( nlc == WHITE ) 1 1 ( nlc != kk ) )continue; 
wpixel( x1,y1,temp_fid ); 

> 

> 

step( &temp_fid,&temp_gid ); 
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} 

} 

} 

printf(*\n\n\t\tFEATURE COUNT IS: %d\fid); 

maplut( IXS, fYS, NCOL, NROW ); 



step( u,v ) /* counter step routine *1 

int *u, *v ; 

{ 

*u +=1; 

if ( *u >= HIGH ) 

{ 

*u = 1; 

*v +=1; 

> 



/* 7 

threshitO 

{ 

int ans, c, option; 
unsigned int val; 

while(l) 

{ 

inrtializeO ; 

geti(“\n\n\tENTER NEW THRESHOLD LEVEL (0255) :',&v al); 

setlut(INPLTT,0); 

threshold (val,HIGHEST); 

setlut(RED.O); 

threshold (val, HIGHEST); 

setlut(GREEN.O); 

threshold (val, HIGHEST); 

setlut(BLUE.O); 

threshold (val, HIGHEST); 

setlut(INPUT,0); 

printf(*\n\nSATISFIED WITH THRESHOLD AT %d ?...Yes (y) / No (n)../,val); 
ans = getch(); 

if(ans == T || ans == ’y) break; 

} 

maplut(IXS,IYS,NCOL,NROW) ; 
inrtializeQ; 



!* V 

geti(name.iptr) /* Get integer keyboard input */ 

char *name; 
int *iptr; 

{ 

printff %s \name); 
sc anf( B %d\iptr); 

} 

getf(rname,riptr) f* Get floating point keyboard input */ 

char *rname; 
int *riptr; 

{ 

printff%s Vname); 
scanty %f\riptr); 
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} 

/* V 

startitO /* initial ITEXPC Board setup */ 

{ 



sethdw (REGBASE, M EM BASE, I FLAG) ; 
aclear(IXS p IYS,NCOL,NROW,1 50) ; 
initialize 0 1 

system Cels'); /* Clears monitor screen V 

} 

/* V 

readitQ 

{ 

int errval; 
int ch; 

while(l) 

{ 

printfC\n\n\tENTER IMAGE FILENAME as (DEV:name.lMG)->"); 

fflush(stdin); 

gets(filename.s); 

filename.l = (char)strlen (filenames); 

errval = readim(IXS,IYS, NCOL.NROW, filename, &comline); 

rf( errval == 0 ){ 

printfC\nFILE-> %s \n\nCOMMENT: %s ",filename.s,comline.s); 
break; 

} 

if( errval 1=0) 

{ 

printf("\n\t\aError reading File!!"); 
rf( errval == -2 ) printf ("\n\tFile Not Found!"); 
rf( errval = = -3) printfC\n\t\tBad File Format!"); 
printf("\n\tTry Again?? ...Yes (y) / No (n) "); 
ch = getch(); 

rf( ch == T || ch == ’y’ ) continue; 
rf( ch == ’N’ || ch == ’n’ ) exit(O); 

> 

> 

> 

/* */ 

saveKO 

{ 

int comflag, errval; 
int ch; 

initial izeO; 

while(l) 

{ 

printff\n\n\tENTER FILENAME... ( DEV: name.lMG )->"); 

fflush(stdin); 

gets (filename. s); 

filename.l = (char)strlen(filename.s); 

getiC\n\nEnter File Compression value ?...( 0/1 ) ",&comflag); 

printfC\n\nEnter Image Comment or Return... to continue.. (Max 200 CHAR):..."); 

fflush(stdin); 

gets(comline.s); 

comline.l = (char) strlen(comline.s); 
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If (comflag == 1) printf("\n\n\tStoring Image using Compression - please wait! 
If (comflag == 0) printff\n\n\tStoring Image - plese wait!! "); 

errval = saveimf IXS,IYS,NCOL,NROW l comflag l filename,comline ); 



if( errval == 0) 

{ 

printf("\n\nlmage Save completed SAT.."); 
break; 

> 

if( errval != 0 ) 

{ 

piintf("\n\aError saving File!!"); 

If( errval == -1 ) printf( "^Insufficient Disk space" ); 
printf("\n\tTry Again?? ...Yes (y) / No (n) "); 
ch = getchO; 

rf( ch == V 1 1 ch == y ) continue; 
if( ch == ’N’ 1 1 ch == 'n' ) break; 

} 

} 

} 

/* */ 
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APPENDIX G 



PROGRAM FILE: sizeit.c 



r 

TITLE : Image FEATURE SIZING Algorithm 
FILENAME : sizeit.c 



LAST 

MODIFIED : 11/26/1988 



PURPOSE : 

Uses an existing or saved image that has been 
thresholded and then ID’d with the ID Algorithm. The output 
of this program is a tabular output of the calculated Area, 
X-chord, and Y-chord that is suitable for input to a 
Statistical Analysis Program. Output file is written to 
default drive as SIZE.DAT. This file should be saved under 
different name each time created as subsequent program 
execution will overwrite the file without warning. 

( Background == WHITE ; Feature == BLACK ) 

V 



#include "thesis.h" 

/* Initial Variable Setup/Definitions */ 

LS20 filename; /* LS20 & LS200 are required for ItexPC */ 

LS200 comline; /* readimO and saveimO library routines */ 

mainO 

{ 

int flag; 

printff\n\aUSE IMAGE processed by the FEATURE ID program!!"); 

printff\n\tReady to Load ID’d image from disk?. ..Yes (y) / No (n) "); 
flag=getch(); 

if (flag == T|| flag == VO 

{ 

startrtO ; 
re ad it 0; 

} 

if( flag == ’N’ 1 1 flag == ’n’ ) exit(O); 

in'rtializeO; 

sizeJtO; 

printff\n\tSave image to Disk File?. ..Yes (y) / No (n) "); 
flag=getchO; 

if (flag == T 1 1 flag == ’y*) savertO; 

} 

I* V 
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size_itO 

{ 

register Int nn, x, y, xa, ya, xb, yb; 

int j, ip, ipa, ipb, xl, yl, x2, y2, maxfl, num; 
int xmax, ymax, xmin, ymin, maxarea, minarea; 

int xflag=0; 
int yflag=0; 
int aflag=0; 

float area, xchord, ychord; 

float fxmax, fymax, fxmin, fymin, fmaxarea, fminarea; 
float cf_1, cf_2, mag; 

long *aptr; 
long *xptr; 
long *yptr; 
long *stop; 

FILE * featdata; 

geti("\nEnter the number of Features', &num); 
getf("\n Enter Magnification Factor - , &mag); 
geti("\nEnter the Maximum Feature Length', & maxfl); 

/* Set initial test min/max values */ 

xmax = 0; 
ymax = 0; 
maxarea = 0; 
xmin = 512; 
ymin = 512; 
minarea = 512; 

/* Dynamically allocate memory space for the number of 
features entered by the operator. 

Note: This eliminates the need to declare large arrays 

for these values and allows maximum number of features */ 

xptr = ( long *) calloc( num .sizeof(long)); 
yptr = ( long *) calloc( num .sizeof(long)); 
aptr = ( long *) calloc( num .sizeof(long)); 

/* Check for successful! memory allocation */ 

rf( I xptr 1 1 lyptr 1 1 (aptr ) 

{ 

printffViOut of Memory!!! \n'); 

printf(^nYou may have entered too many featuresAn') ; 

printffAnCorrect problem and try again."); 

exit(0); 

> 

I* Calculate Magnification Constants */ 

cf^l = ( SCALE_F ACTOR / mag ); 
cf_2 = ( cfj * cfj ); 

/* Above Conversion factor converts size to microns 
and is related to the input magnification for the 
hologram. The defined " SCALE_FACTOR * constant 
value must be set for the proper value based on the 
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equipment set up during actual Image acquisition. 

The constant term SCALE FACTOR is defined In the 
Header File THESIS.H •/ 

I* Begin Sizing Routine *1 

printff\n\n\tStarting to SIZEI...*); 

nn = 0; /* Feature number counter */ 

for( y = 0 ; y < NROW ; y+ + ) 

{ 

if( nn == num ) break; f* Quit when all features sized */ 
for( x = 0 ; x < NCOL ; x++ ) 

{ 

ip = rpixel(x,y); 

if( (ip == WHITE) 1 1 (ip == BLACK) ) continue; 

iff nn == num )break; J* Quit when all features sized */ 

xl = x - maxfl; f* Set up sizing Box V 

yi = y; 

x2 = x + maxfl; 
y2 = y + maxfl; 

iff xl < LOW) xl = LOW; 
iff yi < LOW) yi = LOW; 
iff x2 > NCOL ) x2 = NCOL; 
iff y2 > NROW ) y2 = NROW; 

forf ya = yi ; ya < y2 ; ya+ + ) 

{ 

for( xa = xl ; xa < x2 ; xa+ + ) 

{ 

ipa = rpixelfxa.ya); 
iff ipa != ip )continue; 
aflag + + ; 
xflag+ + ; 

} 

iff xflag > *(xptr + nn ) ) *(xptr + nn) = xflag; 
xflag = 0; 

} 

*(aptr + nn) = aflag; 
aflag = 0; 

forf xb = xl ; xb < x2 ; xb++) 

{ 

forf yb = yi ; yb < y2 ; yb++ ) 

{ 

ipb = rpixelfxb.yb); 
iff ipb 1= ip ) continue; 
yflag+ + ; 

wpixelfxb.yb, BLACK); 

> 

iff yflag > *(yptr + nn ) ) *(yptr + nn) = yflag; 
yflag = 0; 



/* Calculate Min/Max values */ 

xmax = maxf xmax, *(xptr + nn) ); 
ymax = maxf ymax, *(yptr + nn) ); 
maxarea = maxf maxarea, *(aptr + nn) ) ; 
xmin = minf xmin, *(xptr + nn) ) ; 
ymin = minf ymin, *(yptr + nn) ) ; 
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minarea = min( minarea, *(aptr + nn) ) ; 



nn + + ; f* Increment counter *f 

printf("\n\n\tFeature %d complete.. .\nn); 

} 

} /*.... end initial x/y for loop ....*/ 

/* output data section */ 

printf("\n\a Sending output data to screen and FILE SIZE.DAT !"); 

I* waitO; V 

rf( ( featdata = fopenfsize.daf.V) ) == NULL ) 

{ 

printffCANNOT Open output file %s\n\ *featdata); 
return(l); 

} 

printfC ID NO AREA X-Width Y-Width \n"); 

for( j = 0 ; j < num ; j++ ) 

{ 

fprintf( featdata,*%10d %10.3f %10.3f %10.3f\n*. j+1, *(aptr + j) * cf_2, 

*(xptr + j) * cf_1, *(yptr + j) * cf_1 ); 

printf(An%10d %10.3f %10.3f %10.3f \ j+1, *(aptr + j) * cf_2, 

*(xptr + j) * cf 1 , *(yptr + j) * cf_1 ); 

} 

I* This section can be used to put total information on bottom of data file 
if desired. File written now to be used with STATGRAPHICS */ 

r fprintf (featdata An\n\nMax X-chord= %f Max Y-chord= %f Max Area= %f \ 

xmax*cf_1 ,ymax*cf_1 ,maxarea*cf_2); 

fprintf(featdata,"\nMin X-chord= %f Min Y-chord= %f Min Area= %f 
xmin*cf_1 ,ymin*cf_1 , minarea* cf_2); */ 

printf("\n\n\nMax X-chord= %f Max Y-chord= %f Max Area= %f ", 
xmax*cf_1 ,ymax*cf_1 ,maxarea*cf_2); 

prirrtffViMin X-chord= %f Min Y-chord= %f Min Area= %f ", 
xmin*cf_1 ,ymin*cM ,minarea*cf_2); 

/* Close open file and free allocated memory */ 

fclose(featdata); 

free(xptr); 

free(yptr); 

free(aptr); 



f* V 

geti(name,iptr) /* Read/Get integer keyboard input V 

char *name; 
int *iptr; 

{ 

printf(" %s :",name); 
scanf(" %d",iptr); 

} 

getf(rname,riptr) f* Read/Get real (float) keyboard input */ 
char *rname; 
int *riptr; 

{ 
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} 



printff %s : (Float) ',mam©); 
scant (* %f,riptr); 



/* V 

wa'rtO /* Pause routine used to prompt operator */ 

{ 

printf("\n\n Press Return to continue"); 

fflush(stdin); 

getcharQ; 



> 

I* */ 

startrtO /* initial ITEXPC Board setup */ 

{ 



sethdw(REGBASE,MEMBASE,IFLAG) ; 
aclearflXS.IYS.NCOUNROW.1 50) ; 
initializeO; 

system("cls"); /* Clears monitor screen */ 

> 

/* v 

readitO 

{ 

int errval; 
int ch; 

while(l) 

{ 

pr i ntf f\n\n\tE NTE R IMAGE FILENAME as (DEV:name.lMG)->"); 

fflush(stdin); 

gets(filename.s) ; 

filename. I = (char) strlen (filenames); 

errval = readim((XS,rY'S,NCOL l NROW,filename 1 &comline); 

rf( errval == 0 ){ 

printf("\nFILE-“> %s \n\nCOMMENT: %s ".filename.s.comline.s); 
break; 

> 

if( errval 1=0) 

{ 

printf(^n\aError reading File! I") ; 
rf( errval == -2 ) printff\n\tFile Not Found!"); 
rf( errval == -3) printff\n\t\tBad File Format!"); 
printff\n\tTry Again?? ...Yes (y) / No (n) '); 
ch = getchO; 

if( ch == T 1 1 ch == Y ) continue; 
rf( ch == "N" 1 1 ch = = ’n’ ) exit(O) ; 

> 

} 

> 

/* */ 

saveitO 

{ 

int comflag, errval; 
int ch; 

initializeO; 

while(l) 
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printfC\n\n\tENTER FILENAME. ..( DEV: name.lMG )->*); 

fflush(stdin); 

gets(filename.s); 

filename.! = (char) strlen (filenames); 

getif\nEnter File Compression value ?...( 0/1 ) \&comflag); 

printff\n\n Enter Image Comment or Return ..to continue.. (Max 200 CHAR);...'); 

fflush(stdin); 

gets(comline.s); 

comline.l = (char) strlen(comline.s); 

H (comflag == 1) printff\n\tStoring Image using Compression - please wart!!'); 
if (comflag == 0) printffVAtStoring Image - plese waitl! "); 

errval = saveim( DCS, IYS,NCOL,NROW, comflag, filename, comline ); 



if ( errval = = 0) 

{ 

printff\n\tlmage Save completed SAT..*); 
break; 

> 

if ( errval != 0 ) 

{ 

printff\n\aError saving File!!*); 

if( errval == -1 ) printff "\n Insufficient Disk space' ); 

printff\n\tTry Again?? ...Yes (y) / No (n) "); 

ch = getchO; 

if( ch == *Y* 1 1 ch = = y ) continue; 
jf( ch == ’N’ 1 1 ch == ’n* ) break; 

} 

} 

> 

/* V 
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APPENDIX H 

PROGRAM FILE: Istat.c 



/* 

TrTLE : Local Statistical Filter Algorithm 

FILENAME : lstat.c 



LAST 

MODIFIED : 12/9/1988 



PURPOSE : 

Provides for filtering an image using the local statistical 
algorithm. This program processes the image by using a local 
5x5 array of pixels to calculate a statistical value for the 
central pixel of the local array. 

The area of filtering for the image is controlled by the 
operator by keyboard input of number of ROWS and COLS to 
process. Program will allow for a Maximum input of 480 rows 
and 512 cols with a minimum of 1 row and 1 column. Program 
also requests the image Std. Dev. as calculated by Image Action 
Program. The resulting filtered image pixel value is stored on 
disk in a virtual array, that is written to the image processing 
screen when all calculations are complete. DISK SPACE REQUIRED 
for virtual array storage is 250KB minimum available on the disk 
at start of the routine. 

Total calculation time for a 512x480 image is 15.5 min. 

V 



#include “thesis. h' 

LS20 filename; /* LS20 & LS200 are required for ItexPC */ 
LS200 comline; /* readimO and saveimO library routines */ 

int row, col, srow, scol, times; 
float dev; 

mainO 

{ 

int flag; 

extern int row, col, times; 
extern float dev; 

printf(“\n\tReady to Load IMAGE FILE?... Yes (y) / No (n) “); 
flag=getch(); 

if (flag == T|| flag == V*) 

{ 

startitO; 

readitO; 

} 
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getif\n\nEnter number of rows (Max=460) to use In SPECKLE CALC : *,&srow); 
getif\nEnter number of columns (Max=512) to us In SPECKLE CALC : \&scol); 
getif\n\n\nEnter number of filter Iterations to do...\&time$); 
getif\nEnter the number of IMAGE rows (MAX = 460) to FILTER : \&row); 
getifViEnter the number of IMAGE cols (MAX = 512) to FILTER : ',&col); 
getf("\n\nEnter the standard deviation for the Image: \&dev); 

do 

{ 

speckleO; 

Istato ; 

} 

while( -times > 0 ); 

speckleO; 

printffVnSave Image to Disk File?..Yes (y) / No (n) "); 
fiag=getchO; 

If (flag == T\ | flag == V*) saveltQ; 



} 

I* V 

IstatO 

{ 

VACB *item_array; 
long x,y; 

extern int row, col, times; 
extern float dev; 

register int m, n, nn; 

int pixel; 

float suml , 6um2, svar, Ivar, tvar; 
float Imean, Imean2, stddev, stddev2, k; 

int ldata[SNUM]; /* SNUM defined in Header file */ 

I* create a virtual array for the array of filtered 

pixel values the size of the image ( 512x460 ) */ 

in'ft_ v _array CITEMS.VAR'.sizeoffltems) .NULL) ; 

f* open the virtual array, reserve buffer space for 100 elements */ 

item_array = open_v_array(*ITEMS.VAR\100); 

/* start local statistic fitter routine */ 

if( row > NROW ) row = NROW; 
lf( col > NCOL ) col = NCOL; 

stddev = dev / 255.0 ; f* Normalize deviation value V 

stddev2 = stddev * stddev ; 

printff\n\tFilter Running -> %d Runs to go!!...\times-1); 

for ( y = 2 ; y < row - 2 ; y + +) 

{ 

for ( x = 2 ; x < col - 2 ; x++) 
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{ 

pixel = rpixel( x,y ); 
suml = 0; 
sum2 = 0; 
nn = 0; 

for( n=y-2;n<y + 3;n++) 

{ 

for( m = x- 2;m<x + 3; m++ ) 

{ 

ldata[nn] = rpixel( m,n ); 
suml + = ldata[nn]; 
nn++; 

} 

} 

Imean = suml / SNUM ; 

Imean2 = Imean * Imean ; 

for( nn = 0 ; nn < SNUM ; nn+ + ) 

{ 

svar = (ldata[nn] - Imean ) * ( ldata[nn] - Imean ); 
sum2 + = svar; 

> 

Ivar = sum2 / SNUM ; 

tvar = fabs( ( (Ivar + Imean2) / (stddev2 + 1) ) - Imean2 ) 
k = tvar / ( ( stddev2 * Imean2 ) + tvar ) ; 

G( x,y ) = (int) ( Imean + k * ( pixel - Imean ) ) ; 

> 

> 

/* write contents of the G array to image screen *1 

printffVAaXtWriting filtered image to screen..."); 

for( y = 2 ; y < row - 2 ; y++ ) 

{ 

for( x = 2 ; x < col - 2 ; x+ + ) 
wpixel( x,y.G(x.y) ); 

} 

I* close the virtual array V 
close_v_array (item_array) ; 

> 

I* V 

geti(name.iptr) I* Get Integer keyboard input */ 

char *name; 
int *iptr; 

{ 

printff %s '.name); 
scanf(" %d",iptr); 

} 

getf(rname,riptr) I* Get floating point keyboard input */ 

char *rname; 
int *riptr; 

{ 

printf(*%s ".rname); 
scanff %f,riptr); 

> 

/* */ 
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f* Initial fTEXPC Board setup */ 



start ItO 

{ 

sethdw(REGBASE,MEMBASE,IFLAG) ; 
aclearfIXS.fYS.NCOL.NROW.I 50) ; 
initial izeO; 

systemCcIs"); /* Clears monitor screen */ 

} 

/* V 

read it 0 

{ 

int errval; 
int ch; 

while(l) 

{ 

printfC\n\n ENTER IMAGE FILENAME as (DEV:name.lMG)->*); 

fflush(stdin); 

gets(filename.s); 

filename.l = (char)strlen(filename.s); 

errval = readimfIXS.IYS.NCOL.NROW.filename.&comline); 

if( errval == 0 ){ 

printf(*\nFILE-> %s \n\nCOMMENT: %s '.filename.s.comline.s); 
break; 

> 

if( errval != 0 ) 

{ 

printffVAaError reading File! !") ; 
if( errval == -2 ) printff\n\tFile Not Found!*); 
if( errval == -3) printf("\n\t\tBad File Format!*); 
printffTry Again?? ...Yes (y) / No (n) "); 
ch = getchO; 

if( ch == T 1 1 ch == y ) continue; 
if( ch == 'N' 1 1 ch == 'n 1 ) exit(O); 

> 

> 

> 

/* */ 

saveitO 

{ 

int comflag .errval; 
int ch; 

initial izeO; 

while(l) 

{ 

printf(*\n\nENTER FILENAME... ( DEV; name.lMG )->*); 

fflush(stdin); 

gets(filename.s); 

filename.l = (char)strlen(filename.s); 

getif^nEnter File Compression value ?...( 0/1 ) ".Acomflag); 

printf(*\n\n Enter Image Comment or Return ..to continue.. (Max 200 CHAR);...*); 

fflush(stdin); 

gets(comline.s); 

comline.l = (char) strlen(comline.s); 

if (comflag == 1) printff\n\nStoring Image using Compression - please wait!!*); 
H (comflag == 0) printf(*\nStoring Image - plese wait!! *); 
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errval * *aveim( (XS,IYS t NCOL, NROW, comflag, filename, comline ); 

rf( errval = = 0) 

{ 

printf(*\nlmage Save completed SAT.."); 
break; 

> 

i<( errval 1=0) 

{ 

printf(*\n\aError saving Filell*); 

vf( errval -= -1 ) printf( "\nlnsufficient Disk space* ); 

printffVAtTry Again?? ...Yes (y) / No (n) *); 

ch = getchO; 

if( ch == T 1 1 ch == V* ) continue; 
if( ch == 'N* 1 1 ch == ’n’ ) break; 

} 

} 

} 

I* V 

speckleO 

{ 

extern int srow, scol; 
int m, n, x, y, nn; 
long smax, smin; 
long sdata[NUM]; 

unsigned long deviation, ssum; 

float smean, slocal, stotal, tot; 
float spklindex; 

K( srow >= NROW ) srow = NROW; 
rf( scol >= NCOL ) scol = NCOL; 

tot = (long)srow * scol; 

I* Commence calculation */ 

printf("\n\n\tCalculating Speckle INDEX...*); 

for ( n = 1 ; n < srow-2 ; n++) 

{ 

for ( m = 1 ; m < scol-2 ; m++) 

{ 

smax = 0; 
smin = 260; 
ssum = 0; 
nn = 0; 

for( y = ( n-1 ) ; y < (n+2) ; y++ ) 

{ 

for( x = ( m-1 ) ; x < ( m+2 ) ; x++) 

{ 

sdata[nn] = rpixel(x,y); 

ssum += (long)sdata[nn]; 

rf( smax < sdata[nn] ) smax = sdata[nn]; 

if( smin > sdata[nn] ) smin = sdata[nn]; 

nn+ + ; 

> 

> 

deviation = smax - smin; 
rf( ssum == 0 ) smean = 1; 
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smean = ssum / NUM ; 
slocal = (float) deviation / smean; 
stotal += slocal; 

> 

} 

spklindex = stotal / tot; 

printff\n\n\tThe Calculated %d by %d speckle index is %f # ,srow,scol, spklindex) 



/* V 



I* Virtual Array Access Routines */ 



int init_v_array (filename, rec_size, filchar) 
char ^filename, filchar; 
int rec_size; 

{ 

long size; 

FILE *f; 

f = fopen(filename,*wb'); 
if (f != NULL) 

{ 

size = 0; 
fwrite(&size,4,1,f); 
fwrite(&rec_size,2,1 ,f) ; 
fwrite(&filchar,1,1,f); 
fclose(f); 
return(l); 

} 

else 

return (NULL); 

} 



f* write array size of zero */ 



r 

I* 

I* 



r and array element size 
and fill char */ 

to file header */ 



VACB *open_v_array (filename, buffer_size) 
char ^filename; 
int buffer size; 

{ 

VACB *v_array; 
char A buf_ptr; 
int i; 

char filchar; 



l* allocate virtual array control block */ 

v array = (VACB *) malloc(sizeof(VACB)); 
if (v_array == NULL) return(NULL); 

/* open virtual array file */ 

v array- >file = fopen(filename,*r+b"); 
if" (v_array->file == NULL) 

{ 

free (v_ar ray) ; 
return (NULL); 

}; 



r get array size and element size for control block */ 

fread(&v_array-> size, 4,1 ,v_array->file) ; 
fread(&v_array->elsize,2,1,v - array->file); 
f read(&filchar, 1 , 1 ,v_array- >file) ; 
v_array->buf_elsize = v_array->elsize + 4; 

/* allocate buffer *1 
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v_arr ay- > buffer = (char *) malloc(v_array->buf_elslze * (buffer_size + 1)); 
if (v_array-> buffer == NULL) { 
fclo8e(v_array->file) ; 
free(v_array); 
return (NULL); 

}; 

v_array->buf_size = buffer_size; 

/* set up blank_rec using the fill character in array header */ 
f* for initializing new array elements V 

buf_ptr = v_array-> buffer + v_array->buf_elsize * v_array->buf_size; 

v_array->blank_rec = bufptr + 4 ; 

for (i = 0; i < v_array->buf_elsize; i++) 

*buf_ptr+ + = filchar; 

/* set record index negative for all buffer elements */ 

buf_ptr = v_array-> buffer; 

for (i = 0; i < v_array->buf_$ize; i++) 

{ 

*((long *)buf_ptr) = -1L; 

buf ptr + = v_array->buf_elsize; 

}; 

return (v_array); 



void close_v_array(v_array) 

VACB *v array; 

{ 

int i; 

char *buf_ptr; 

long recindex, file_offset; 

buf_ptr = v_array-> buffer; 

I* flush buffer */ 

for (i = 0; i < v_array->buf_size; i++) 

/* check each element index */ 

I* if element present; write it to disk */ 

rec index = *((long *)buf_ptr); 
if (recjndex >= 0) 

{ 

file_offset = header + recjndex * v_array->elsize; 

fseek (v_array->file,file_offset,0) ; 

fwr’rte(buf ptr + 4, v array- >elsize,1, v_ar ray- > file); 

}; 

buf_ptr + = v_array->buf_elsize; 



void ‘access j/jec(v_array, index) 
VACB *v_array; 
long index; 

{ 

char *buf_ptr; 
int bufjndex; 



} 



}; 

free(v_array- > buffer) ; 
fclose(v_array->file) ; 
free(v_array); 



f* de-allocate buffer */ 
f* close array file */ 

I* de-allocate VACB •/ 



} 
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long recjndex, tempjndex; 



/* calculate buffer address of referenced element */ 

bufjndex = Index % v_array->buf_size; 

buf_ptr = v_array-> buffer + bufjndex * v_array->buf_elsize; 

recjndex = *(long *)buf_ptr; 

/* If element present, return its buffer address */ 

If (recjndex - = Index) retum(buf_ptr + 4); 

/* If element doesn’t exist, extend the file V 

If (index >= v_array->size) 

{ 

fseek(v_array-> file, 0,2) ; 

for (tempjndex = v_array->size; tempjndex+ + <= index; ) 
fwrKe(v_airay-> blank jec, v_array->elslze, 1, v_array*>file); 
v_array->size = index + 1; 
fseek(v_array-> file, 0,0) ; 
fwrite(&v_array->size t 4, 1, v_array->file); 

}; 



/* If buffer slot is occupied by another e.ement, */ 

I* save it to disk V 

If (recjndex >= 0) 

{ 

fseek(v_array->file, recjndex * v_array->elsize + header, 0); 
fwrite(buf_ptr + 4, v_array->elsize, 1, v_array->file); 

}; 

I* read referenced element into buffer slot V 

fseek(v_array->file, index * v_array->elsize + header, 0); 
fread(buf_ptr + 4, v_array->elsize, 1, v_array->file); 

*((long *)buf_ptr) = index; 

/* return address of element in buffer */ 

return (buf_ptr + 4); 

} 

/* */ 
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APPENDIX I 



PROGRAM FILE: 2$lgma.c 



/* 

TITLE : 2SIGMA Statistical Filter Algorithm 

FILENAME : 2sigma.c 



LAST 

MODIFIED : 12/9/1988 



PURPOSE : 

Provides for filtering an image using a 2sigma statistical 
algorithm. This program processes the image by using a local 
5x5 array of pixels to calculate a statistical value for the 
central pixel of the local array. The algorithm will eliminate 
pixels from the summation that are greater than 2sigma in 
value from that of the central pixel. 

The area of filtering for the image is controlled by the 
operator by keyboard input of number of ROWS and COLS to 
process. Program will allow for a maximum input of 480 rows 
and 512 cols with a minimum of 1 row and 1 column. Program 
also requests the image standard deviation as obtained using 
Image Action Program. The resulting filtered pixel value is 
stored on disk in a virtual array and written to screen later. 

Total calculation time for a 512x480 image is 5.7 min. 

DISK SPACE REQUIRED for virtual array storage is 250 KB 
minimum available on the disk at start of the routine. 

V 

#include “thesis.h" 

LS20 filename; f* LS20 & LS200 are required for ItexPC *f 
LS200 comline; /* readimO and saveimO library routines */ 

int row, col, srow, scol, times; 
float dev; 

main() 

{ 

int flag; 

extern int row, col, times; 
extern float dev; 

printffVAtReady to Load IMAGE FILE from disk?. ..Yes (y) / No (n) '); 
flag = g etch 0; 

if (flag = = Y 1 j | flag == ’y*) 

{ 

start rtO; 
readftQ; 



} 
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geti(*\n\nEnter number of rows (Max=480) to use in SPECKLE CALC : ',&srow); 

getif\nEnter number of columns (Max=512) to us in SPECKLE CALC : \&scol); 

getiOn\n\nEnter number of filter iterations to do...',&times); 

getifViEnter the number of IMAGE rows (MAX = 480) to FILTER : *,&row); 

getifViEnter the number of IMAGE cols (MAX = 512) to FILTER : \&col); 

getf(*\n\nEnter the standard deviation for the Image: \&dev); 

do 

{ 

speckleO; 

sigmaO; 

while( -times > 0 ); 
speckleO; 

printff\nSave Image to Disk File?..Yes (y) / No (n) "); 
flag =g etch 0; 

if (flag = = V| | flag == V 1 ) aaveitQ; 



} 

/* *f 

sigmaO 

{ 

VACB *rtem_array; 
long x,y; 

extern int row, col, times; 
extern float dev; 

register int m, n, nn; 

int pixel, suml, sum2, delta, high, low; 

float stddev, hvar, Ivar; 

int IdatafSNUM]; /* SNUM defined in Header file *1 

I* create a virtual array for the array of filtered 

pixel values the size of the image ( 512x480 ) */ 

i nit^v^arrayC ITEMS. VAR', sizeof (items) , NU LL) ; 

/* open the virtual array, reserve buffer space for 100 elements V 

item_array = open_v_array(TTEMS.VAR , ,100); 

I* start local statistic filter routine V 

if( row > NROW ) row = NROW; 
if( col > NCOL ) col = NCOL; 

stddev = dev J 255.0 ; /* Normalize deviation value */ 

hvar * 1. + 2. * stddev ; /* find high and low sigma */ 

Ivar = 1.-2.* stddev ; 

printf(*\n\n\tFilter Running -> %d Runs to go!L..\times-1); 

for ( y = 2 ; y < row - 2 ; y++) 

{ 

for ( x = 2 ; x < col - 2 ; x++) 

{ 
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pixel = rpixel ( x ( y ); 
suml = 0; 
sum2 = 0; 
nn = 0; 
delta - 0; 

high = (int)(hvar * pixel); 
low = (int)(lvar * pixel); 

for( n = y- 2;n<y + 3;n++) 

{ 

for( m = x- 2;m<x + 3; m++ ) 

{ 

ldata[nn] = rpixel( m,n ); 

rf( ( low <= Idatafnn] ) && ( ldata[nn] <= high ) ) 

{ 

suml += ldata[nn]; 
delta+ + ; 

} 

nn+ + ; 

} 

} 

if(delta < = 2) /* correct shot noise - 4 neighbor average */ 

{ 

sum2 = ( r pixel (x,y-1 ) + r pixel (x p y + 1 ) + rpixel (x-1 ,y) + rpixel (x+ 1 ,y) ); 

G( x p y ) = (pixel + sum2) / 5; 

continue; 

> 

G( x,y ) = suml / delta ; 

> 

} 

printf("\n\aWriting filtered image to screen...'); 

for( y m 2 ; y < row - 2 ; y++ ) 

{ 

for( x = 2 ; x < col - 2 ; x++ ) 
wpixel( x p y p G(x p y) ); 

} 

/* close the virtual array */ 
close_v_array (item_array) ; 

} 

I* V 

geti(name.iptr) I* Get integer keyboard input */ 

char *name; 
int *iptr; 

{ 

prirrtff %s '.name); 
scanfC %d' p ipti^; 

} 

getf(rname,riptr) f* Get floating point keyboard input */ 

char *rname; 
int *riptr; 

{ 

printfC%s \rname); 
scanff %f,riptr); 

> 
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f* V 

startrtO /* initial ITEXPC Board setup */ 

{ 



sethdw(REG BASE, MEM BASE, I FLAG) ; 
aclear(IXS,IYS,NCOL,NROW,150); 
inrtializeO ; 

system (■els*); /* Clears monitor screen *f 

} 

f* V 

read itO 

{ 

int errval; 
int ch; 

while(l) 

{ 

printf("\n\nENTER IMAGE FILENAME as (DEV:name.lMG)->"); 

fflush(stdin); 

gets(filename.s); 

filename.l = (char)strlen(filename.s); 

errval = readim(IXS t (YS,NCOL,NROW,filename,&comline); 

if ( errval == 0 ){ 

printfC\nFILE-> %s \n\nCOMMENT: %s , ,filename.s,comline.s); 
break; 

} 

if( errval 1=0) 

{ 

printf(“\n\aError reading File!!'); 
if( errval = = -2 ) printffVAtFile Not Found!"); 
rf( errval = = -3) printff\n\t\tBad File Format!'); 
printffTry Again?? ...Yes (y) / No (n) "); 
ch = getchO; 

rf( ch == T 1 1 ch == y ) continue; 
if( ch == 'N’ || ch == ’n’ ) exit(O); 

} 

} 

} 

/* */ 

saveitO 

{ 

int comflag, errval; 
int ch; 

inrtializeO ; 

while(l) 

{ 

printff\n\nENTER FILENAME.. .( DEV: name.lMG )->"); 

fflush(stdin); 

gets(filename.s); 

filename.l = (char) strlen (filenames); 

geti(^nEnter File Compression value ?...( 0/1 ) " ( &comflag); 

printff \n\n Enter Image Comment or Return ..to continue.. (Max 200 CHAR):...’); 

fflush(stdin); 

gets(comline.s); 

comline.l = (char) strlen(comline.s); 

if (comflag == 1) printff\n\nStoring Image using Compression - please wart!!"); 
if (comflag == 0) printff\nStorlng Image - plese wait!! "); 
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errval = saveim( IXS.IYS.NCOUNROW.comflag.filename.comline ); 



if( errval == 0) 

{ 

printff\nlmage Save completed SAT.."); 
break; 

> 

if ( errval 1=0) 

{ 

printff\n\a Error saving File! I"); 

if( errval = = -1 ) printf( *\nlnsufficient Disk space" ); 

printf("\n\tTry Again?? ...Yes (y) / No (n) "); 

ch = getchO; 

lf( ch == V 1 1 ch == y ) continue; 
rf( ch == ’N’ 1 1 ch == 'n‘ ) break; 

} 

} 

> 

/* */ 

speckle)) 

{ 

extern int srow, scol; 
int m, n, x, y, nn; 
long smax, smin; 
long sdata[NUM]; 

unsigned long deviation, ssum; 

float smean, slocal, stotal, tot; 
float spklindex; 

if( srow >= NROW ) srow = NROW; 
if( scol >= NCOL ) scol = NCOL; 

tot = (long)srow * scol; 

/* Commence calculation */ 

printf("\n\n\tCalculating Speckle INDEX...'); 

for ( n = 1 ; n < srow-2 ; n++) 

{ 

for ( m = 1 ; m < scol-2 ; m++) 

{ 

smax = 0; 
smin = 260; 
ssum = 0; 
nn = 0; 

for(y = ( n-1 ) ; y < (n+2) ; y+ + ) 

{ 

for( x = ( m-1 ) ; x < ( m+2 ) ; x++) 

{ 

sdata[nn] = rpixel( x,y ); 

ssum += (long)sdata[nn]; 

if( smax < sdata[nn] ) smax = sdata[nn]; 

if( smin > sdata[nn] ) smin = sdata[nn]; 

nn+ + ; 

} 

} 

deviation = smax - smin; 
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If( Mum == 0 ) smean = 1; 
smean = ssum / NUM ; 
slocal = (float) deviation / smean; 
stotal += slocal; 

} 

} 

spklindex = stotal / tot; 

printf(*\n\n\tThe Calculated %d by %d speckle index is %f \srow,scol, spklindex) 



/* */ 

/* Virtual Array Access Routines */ 



int init_v_array(filename l rec_size l filchar) 
char ‘filename, filchar; 
int rec_size; 



long size; 

FILE M; 

f = fopen(filename l “wb'); 
rf (f 1= NULL) 

{ 

size = 0; 

fwrrte(&size,4,1,f); 

fwrrte(&rec_size,2,1 ,f); 

fwrrte(&filchar,1,1,f); 

fclose(f); 

return (1); 

} 

else 

ret urn (NULL); 

} 



/* write array size of zero */ 

/* and arruy element size */ 
/‘ and fill char */ 

I* to file header */ 



VACB ‘open_v_array (filename, buffer_size) 
char ‘filename; 
int buffer_size; 

{ 

VACB ‘v_array; 
char ‘buf_ptr; 
int i; 

char filchar; 



/* allocate virtual array control block */ 

v_array = (VACB ‘) malloc(sizeof(VACB)); 
rf (v_array == NULL) return(NULL); 

/* open virtual array file */ 

v_array->file = fopen(filename t V+b*); 
if (v_array->file == NULL) 

{ 

free(v_array) ; 
return (NULL); 

}; 



f* get array size and element size for control block ‘/ 

tread (&v_array->size,4,1,v_array->file); 
tread (&v_array->elsize,2,1 ,v_array->file) ; 
tread (&filchar,1 ,1 ,v_array->file) ; 
v_array->buf_elsize = v_array->elsize + 4; 

/‘ allocate buffer */ 
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v_array-> buffer = (char *) malloc(v_array->buf_elsize * (buffer_size + 1)) 
if (v_array-> buffer == NULL) { 
f close (v_array- > f i le) ; 
free(v_array) ; 
return (NULL); 

}; 

v_array->buf_size = buffer_slze; 

/* set up blank_rec using the fill character in array header */ 
r for initializing new array elements *1 

buf_ptr = v_array-> buffer + v_array->buf_elsize * v_array->buf_size; 

v_array->blank_rec = buf_ptr + 4; 

for (i = 0; i < v_array->buf_elsize; i++) 

*buf_ptr++ = filchar; 

I* set record index negative for all buffer elements */ 

buf_ptr = v_array-> buffer; 

for (i = 0; i < v_array->buf_size; I++) 

{ 

*((long *)buf_ptr) = -1L; 
buf_ptr + = v_array->buf_elsize; 

}; 

return(v_array); 



void close_v_array(v_array) 

VACB *v_array; 

{ 

int i; 

char *buf_ptr; 

long recjndex, file_offset; 

buf_ptr = v_array-> buffer; 

I* flush buffer */ 

for (i = 0; i < v_array->buf_size; i++) 

{ 

/* check each element index */ 

/* if element present; write it to disk */ 

rec_index = *((long *)buf_ptr); 
if (recjndex >= 0) 

{ 

file_offset = header + recjndex * v_array->elsize; 

fseek(v_array- >file,file_offset,0) ; 

fwrite(buf_ptr + 4, v_array->elsize,1, v_array->file); 

}; 

buf_ptr + = v_array->buf_elsize; 



void *access_v_rec(v_array, index) 
VACB # v_array; 
long index; 

{ 

char *buf_ptr; 
int bufjndex; 



} 



}; 

free(v_array- > buffer) ; 
fclose(v_array->file) ; 
free(v_array) ; 



I* de-allocate buffer */ 
f* close array file */ 

/* de-allocate VACB */ 



} 
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long recjndex, tempjndex; 



/* calculate buffer address of referenced element */ 

bufjndex = index % v_array->buf_size; 

buf_ptr = v_array-> buffer + bufjndex * v_array->buf_elsize; 

recjndex = *(long *)buf_ptr; 

I* if element present, return its buffer address *1 
rf (recjndex == index) return (buf_ptr + 4); 

f* rf element doesn’t exist, extend the file */ 

if (index >= v_array->size) 

{ 

fseek(v_array->file,0,2) ; 

for (tempjndex — v_array->size; tempjndex+ + <= index; ) 
fwrrte(v_array->blank_rec, v_array->elsize, 1, v_array->file); 
v_array->size = index + 1; 
fseek(v_array->file,0,0); 
fwrite(&v_array->$ize, 4, 1, v_array->file); 

}; 

I* if buffer slot is occupied by another element, */ 

I* save it to disk */ 

if (recjndex >= 0) 

{ 

f$eek(v_array->file, recjndex * v_array->elsize + header, 0); 
fwrite(buf ptr + 4, v array- >elsize, 1, v_array->file); 

}; 

I* read referenced element into buffer slot V 

fseek(v_array->file, index * v_array->el$ize + header, 0); 
fread(buf_ptr + 4, v_array->elsize, 1, v_array->file); 

•((long *)buf_ptr) = index; 

I* return address of element in buffer V 

return (buf_ptr + 4); 

} 

/* */ 
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APPENDIX J 



PROGRAM FILE: geofil.c 



/* 

TITLE : Geometric Filter Algorithm 

FILENAME : geofil.c 



LAST 

MODIFIED : 12/9/1988 



PURPOSE : 

Provides for filtering an image using the geometric filter 
algorithm. This program processes the image and provides 
for the operator to select the number of iterations. 

The area of filtering for the image is controlled by the 
operator by keyboard input of number of ROWS and COLS to 
process. Program will allow for a Maximum input of 480 rows 
and 512 cols with a minimum of 1 row and 1 column. 

The resulting filtered image pixel value is stored on 
disk in a virtual array that is written to the image 
processing screen when all calculations are complete. 

Total calculation time for a 512x480 image is 41.7 min. 

DISK SPACE REQUIRED for virtual array storage is 250 KB minimum 
available on the disk at start of the routine. 

V 



#include thesis. h* 

LS20 filename; /* LS20 & LS200 are required for KexPC */ 

LS200 comline; /* readimO and saveimO library routines */ 

int row, col, srow, scol, times; 

mainO 

{ 

int flag; 

extern int row, col, srow, scol, times; 

printf(“\n\tReady to Load IMAGE FILE?...Ye$ (y) / No (n) ’); 
flag=getch(); 

if (flag == T|| flag == V*) 

{ 

startitO; 

readitO; 

> 

rf( flag == *N* 1 1 flag == 'n’ ) ex'rt(O); 

getif\n\n Enter number of rows (Max=480) to use in SPECKLE CALC : - ,&srow); 
geti(“\nEnter number of columns (Max=512) to us in SPECKLE CALC : *,&scol); 
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getifVArAn Enter number of filter Iterations to do...*,&times); 
getifViEnter the number of IMAGE rows (MAX = 480) to FILTER : \&row); 
getif\nEnter the number of IMAGE cols (MAX = 512) to FILTER : ',&col); 

do 

{ 

speckleO; 

geofilO; 

} 

while( -times > 0 ); 
speckleO; 

printffAnSave image to Disk File?..Yes (y) / No (n) ■); 
flag=getchO; 

if (flag == T| | flag = = V) saveitQ; 



} 

/* V 

geofilO 

{ 

VACB *item_array; 
long x,y; 

extern int row, col, times; 

int fl, f2, gl , g2, g3; 
int pixel, a, b, c, d; 

I* create a virtual array for the array of filtered 

pixel values the size of the image ( 512x480 ) */ 

inrt_v_ arr ay(TrEMS.VAR*,sizeof (items), NULL); 

I* open the virtual array, reserve buffer space for 100 elements */ 

item_array = open_v_array(TnEMS.VAR\100); 

/* start local statistic filter routine */ 

if ( row > NROW ) row = NROW; 
if ( col > NCOL ) col = NCOL; 

for( y = 1 ; y < row ; y++ ) /* Zero the initial array */ 

{ 

for( x = 1 ; x < col ; x++ ) 

{ 

G( x,y ) = NULL; 

> 

> 

printf("\n\n\tFilter Running “> %d Runs after this onel...\n*,times-1); 

a=1; b=0; c=3; d=1; 
printffAnc - 3 w )\ 

tl: for( y = 1 ; y < row-1 ; y++ ) 

{ 

for( x = 1 ; x < coH ; x++ ) 

{ 



fl = rpixel( x,y ); 
f2 = rpixel( x-a,y-b ); 

G( x,y ) = max( fl, min( f2-1, fl+1 ) ); 
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} 



} 

printff step 1 ■); 

for( y = 1 ; y < row-1 ; y++ ) 

{ 

for( x = 1 ; x < coM ; x+ *f ) 

{ 

01 • G(x,y); 

g2 = G(x-a,y-b); 
g3 = G( x+a,y+b ); 

pixel = max( gl, min( g2, min( gl + 1, g3+1 ) ) ); 
wpixel( x,y, pixel); 

} 

} 

printfC step 2*); 

»( d= = 1 ) 

{ 

a= -a; 
b= -b; 
d=0; 
goto tl ; 

} 

rf( d==0 ) 

{ 

<1=1 ; 

goto t2; 

} 

t2: for( y = 1 ; y < row-1 ; y++ ) 

{ 

for( x = 1 ; x < coM ; x++ ) 

{ 

fl = rpixel( x,y ); 

12 = rpixel( x-a,y-b ); 

G( x,y ) = min( fl, max( f2-1, fl +1 ) ); 

} 

> 

printfC step 3‘) ; 

for( y = 1 ; y < row-1 ; y++ ) 

<or( x = 1 ; x < col-1 ; x++ ) 

{ 

91 = G( x,y); 

g2 = G( x-a,y-b ); 

03 = G( x+a,y+b ); 

pixel - max( gl, min( g2, min( gl + 1, g3+1 ) ) ); 
wpixelf x,y, pixel); 

> 

> 

printff step 4‘); 
if(d==1 ) 

{ 

a= -a; 
b= -b; 
d=0; 
goto t2; 

> 

rf( d==0 ) 

{ 

d=1; 
goto t3; 

} 

t3: switch (c) 
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{ 

case 0: printfC\nc = 0 END "); 
break; 

case 3: a=0; 
b=1; 
c=2; 

printffVic = 3 *); 
goto tl ; 

case 2: a=1 ; 
b=1; 
c=i; 

printffVic = 2 "); 
goto tl ; 



case 1: a=1; 

b= -1; 
c=0; 

printf(*\nc = 1 "); 
goto tl ; 

default : break; 

} 

/* close the virtual array */ 
close v_array(item_array); 

} 

r v 

geti(name.iptr) /* Get integer keyboard input */ 

char ‘name; 
int *iptr; 

{ 

printff %s ’.name); 
scanf(* %d*,iptr); 

} 

getf(rname,riptr) /* Get floating point keyboard input */ 

char *rname; 
int *riptr; 

{ 

printff%s , ,rname); 
scanff %f’,riptr); 



} 

r v 

startitO /* initial ITEXPC Board setup */ 

{ 



sethdw(REGBASE,MEMBASE,IFLAG) ; 

aclear(IXS,IYS,NCOL,NROW,150); 

initializeO; 

systemfcls*); /* Clears monitor screen */ 

} 

r */ 

readitO 

int errval; 
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int ch; 



while(l) 

{ 

printf("\n\n\tENTER IMAGE FILENAME as (DEV:name.lMG)->"); 

fflush(stdin); 

gets(filename.s) ; 

filename.l = (char)strlen(filename.s); 

errval = readim(DCS,IYS, NCOL,NROW, filename, &comline); 

if( errval == 0 ){ 

prlntff\nFILE-> %s \n\nCOMMENT: %s ", filenames, comline.s); 
break; 

> 

if( errval 1=0) 

{ 

printf("\n\aError reading Filel I") ; 
ff( errval == -2 ) printff\n\tFile Not Foundl"); 
rf ( errval == -3) printf(Vi\t\tBad File Formatl"); 
printffTry Again?? ...Yes (y) / No (n) "); 
ch = getchO; 

H( ch = = T 1 1 ch == Y ) continue; 

*rf( ch == ’N’ 1 1 ch == ’n* ) exit(O); 

} 

} 

} 

/* */ 

saveitO 

{ 

int comflag, errval; 
int ch; 

inrtializeO; 

while(l) 

{ 

printf("\n\n\tENTER FILENAME... ( DEV: name.lMG )->"); 

fflush(stdin); 

gets(filename.s); 

filename.l = (char)strlen(filename.s); 

geti(^n Enter File Compression value ?...( 0/1 ) \&comflag); 

printff \n\n Enter Image Comment or Return ..to continue.. (Max 200 CHAR):..."); 

fflush(stdin); 

gets(comline.s); 

comline.l = (char) strlen(comline.s); 

if (comflag == 1) printf("\n\nStoring Image using Compression - please wart!!*); 
if (comflag == 0) printf(*\nStoring Image - plese wait!! "); 

errval = saveim( IXS.IYS.NCOUNROW.comflag.filename.comline ); 



if ( errval = = 0) 

{ 

printf(*\nlmage Save completed SAT.."); 
break; 

} 

if ( errval 1=0) 

{ 

printf(*\n\a Error saving File!!"); 

if( errval == -1 ) printf( "\nlnsufficient Disk space" ); 

printf("\n\tTry Again?? ...Yes (y) / No (n) "); 

ch = getchO; 
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if( ch == T 1 1 ch == Y ) continue; 
lf( ch == ’N’ || ch «« 'n’ ) break; 

} 

> 

} 

/* */ 



speckleO 

{ 

extern int srow, scol; 
int m, n, x, y, nn; 
long smax, smin; 
long sdata[NUM]; 

unsigned long deviation, ssum; 

float smean, slocal, stotal, tot; 
float spklindex; 

if( srow >= NROW ) srow = NROW; 
if( scol >= NCOL ) scol = NCOL; 

tot = (long) srow * scol; 

/* Commence calculation */ 

printf("\n\n\tCalculating Speckle INDEX..."); 

for ( n = 1 ; n < srow-2 ; n++) 

{ 

for ( m = 1 ; m < scol-2 ; m++) 

{ 

smax = 0; 
smin = 260; 
ssum = 0; 
nn = 0; 

for( y = ( n-1 ) ; y < (n+2) ; y+ + ) 

{ 

for( x = ( m-1 ) ; x < ( m+2 ) ; x++) 

sdata[nn] = rpixel( x,y ); 

ssum += (long)sdata[nn]; 

if( smax < sdata[nn] ) smax = sdata[nn]; 

if( smin > sdatafnn] ) smin = sdata[nn]; 

nn+ + ; 

} 

} 

deviation = smax - smin; 
rf( ssum == 0 ) smean = 1; 
smean = ssum / NUM ; 
slocal = (float) deviation / smean; 
stotal += slocal; 

} 

} 

spklindex = stotal / tot; 

printfp\n\n\tThe Calculated %d by %d speckle index is %f \srow,scol, spklindex); 



/* V 

I* Virtual Array Access Routines */ 

int inlt_v_array(filename,rec_8ize,fiichar) 
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char ^filename, filchar; 
int rec_size; 

{ 

long size; 

FILE *f; 

f = fopen (filename, *wb"); 

If (f != NULL) 

size = 0; 
fwrite(&size,4,1,f); 
fwrite(&rec_slze,2,1 ,f); 
fwrrte(&filchar,1 ,1 ,f); 
fclose(f); 
return(l); 

> 

else 

return (NULL); 

} 

VACB *open_v_array (filename, buffer_size) 
char ^filename; 
int buffer_size; 

{ 

VACB *v_array; 
char *buf_ptr; 
int i; 

char filchar; 

/* allocate virtual array control block */ 

v_array = (VACB *) malloc(sizeof(VACB)); 
if (v_array == NULL) return(NULL); 

/* open virtual array file */ 

v_array->file = fopen (filename, 'r+b'); 
if (v_array->file == NULL) 

{ 

free(v_array) ; 
return (NULL); 

}; 



/* write array size of zero */ 
f* and array element size */ 
f* and fill char */ 

I* to file header */ 



/* get array size and element size for control block */ 

fread(&v_array->size,4,1,v_array->file); 
freed (&v_array-> elsize, 2,1 ,v_array->file) ; 
fread(&filchar,1,1,v_array->file); 
v_array->buf_elsize = v_array-> elsize + 4; 

/* allocate buffer *1 

v_array-> buffer = (char *) malloc(v_array->buf_elsize * (buffer_size + 1)); 
if (v_array-> buffer == NULL) { 
fclose(v_array->file) ; 
free(v_array) ; 
return (NULL); 

}; 

v_array->buf_size = buffer_size; 

/* set up blank_rec using the fill character in array header */ 

/* for initializing new array elements */ 

buf_ptr = v_array-> buffer + v_ar ray- >buf_el size * v_array->buf_size; 
v_array->blank_rec = buf_ptr + 4; 
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for (i = 0; I < v_array->buf_elsize; I + +) 

*buf_ptr++ = filchar; 

/* set record index negative for all buffer elements */ 

buf_ptr = v_array-> buffer; 

for (i = 0; I < v array- >buf_size; I++) 

{ 

*((long *)buf_ptr) = -1L; 
buf_ptr + = v_array->buf_elsize; 

}; 

return (v_arr ay); 



void close_v_array(v_array) 
VACB *v_array; 

{ 

Int i; 

char *buf_ptr; 

long recjndex, file_offset; 

buf_ptr = v_array-> buffer; 

/* flush buffer */ 



for (i = 0; i < v_array->buf_size; i++) 

{ 

/* check each element index */ 

/* if element present; write it to disk */ 

recjndex = *((long *)buf_ptr); 
if (rec_index >= 0) 

{ 

file_offset = header + recjndex * v_array->elsize; 

fseek (v_array- > file,file_offset,0) ; 

fwrite(buf_ptr + 4, v_array->elsize,1 , v_array«>file); 

>; 

buf_ptr += v_array->buf_elsize; 

>; 

free(v_array-> buffer); /* de-allocate buffer */ 

fclose(v_array->file); /* close array file */ 

free(v_array) ; /* deallocate VACB *1 



void *access_v_rec(v_ar ray, index) 
VACB *v_array; 
long index; 

{ 



char *buf_ptr; 
int buf_index; 

long recjndex, tempjndex; 



r calculate buffer address of referenced element */ 

bufjndex = Index % v_array->buf_size; 

buf_ptr = v_array-> buffer + bufjndex * v_array->buf_elsize; 

recjndex = *(long *)buf_ptr; 

/* if element present, return its buffer address */ 
if (recjndex == Index) return(buf_ptr + 4); 

/* if element doesn’t exist, extend the file */ 
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if (index >= v_array->slze) 

{ 

teeek(v_array->file,0,2) ; 

for (tempjndex = v_array->size; temp_index+ + <= index; ) 
v fwrite(v_array->blank_rec ( v_array->elsize, 1, v_array->file); 

v_array->size = index + 1; 
faeek(v_array->file,0,0) ; 
fwrrte(&v_array->stze, 4, 1, v_array->file); 

>; 

/* if buffer alot is occupied by another element, */ 

/* save it to disk */ 

if (recjndex >= 0) 

{ 

fseek(v_array->file, recjndex * v_array->elsize + header, 0); 
fwrite(buf_ptr + 4, v array- >elsize, 1, v_array->file); 

>; 

I* read referenced element into buffer slot */ 

fseek(v_array->file, index * v_array->elsize + header, 0); 
fread(buf_ptr + 4, v_array->elsize, 1, v_array->file); 

*((long *)buf_ptr) = index; 

/* return address of element in buffer V 

return(buf_ptr + 4); 

> 

i* */ 
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APPENDIX K 



FORTRAN AND C COMPILER OPERATIONS 

The following DOS batch files list the FORTRAN and C compiler directives used to compile 
all the FORTRAN and C programs using the optimizing compilers. All of these batch files use 
the recommended command line compile directives found in the ITEX/PC programmers manual 
[Ref. 6]. 

The FORTRAN batch files consist of FCOMP.BAT, FLINK.BAT, FPACK.BAT, and FORALL.BAT. 
The programs should be located in a directory named FORTRAN along with the program to 
be compiled. The batch files are executed using a DOS command that includes the name of 
the file to be compiled or linked. An example of its use is as follows: 

FCOMP newid 

which will compile the file newid.for located in the FORTRAN directory on the DOS path. 

The other batch files operate in a similar manner. The file FORALLBAT combines the compile, 
link, and pack files into one evolution. This file is useful for compiling and linking programs 
known to be free of errors with just one DOS command. 

File FCOMP.BAT 

PATH=\fortran;c:\fort4\src;c:\fort4;c:\fort4\bin;c:\fort4\lib;c:\fort4\include; 

SET INCLUDE=c:\fort4\include 
SET LIB=c:\F0RT4\LIB 
SET TMP=c:\fort4\temp 

FL /c /Fsd:\fortran\%1 /Fcd:\fortran\%1 %1.for 



File FLINK.BAT 

PATH=\fortran;c:\fort4\src;c:\fort4;c:\fort4\bin;c:\fort4\lib;c:\fort4\include; 
SET INCLUDE=c:\fort4\include 
SET LIB=c:\FORT4\LIB 
SET TMP= c:\fort4\temp 

link %1 %2 %3,„rrEXPC.uB fortran.lib; 



File FPACK.BAT 

PATH=\fortran;c:\fort4\src;c:\fort4;c:\fort4\bin;c:\fort4\lib;c:\fort4\include; 

SET INCLUDE=c:\fort4\include 

SET LIB=c:\F0RT4\LIB 

SET TMP=c:\fort4\temp 

exepack %1.exe %1.tmp 
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COPY %1.TMP %1 .EXE 
ERASE %1.TMP 

File FORALL.BAT 

PATH=\fortran;c:\fort4\src;c:\fort4;c:\fort4\bin;c:\fort4\lib;c:\fort4\include;SET 

INCLUDE=c:\fort4\include 

SET LIB=c:\FORT4\UB 

SET TMP=c:\fort4\temp 

FL /c /Fsd:\fortran\%1 /Fcd:\fortran\%1 %1.for 

LINK %1 %2 %3,„rrEXPC.LIB FORTRAN. LIB; 

EXEPACK %1.EXE %1.TMP 
COPY %1.TMP %1.EXE 
ERASE %1.TMP 
ERASE *.COD 
ERASE *.LST 
ERASE *.OBJ 
ERASE *.MAP 
ERASE *.BAK 



The C compiler file ACOMP.BAT performs essentially the same evolutions as the FORTRAN 
files but only one file is necessary to accomplish all of the evolutions. Use of this file requires 
the program to be compiled to be in the same directory as the ACOMP.BAT file and to be on 
the path specified. The program will be compiled and, if successful, the link evolution will also 
occur. The defaults used with the compile command line achieve the generation of packed 
code automatically. All C programs in this thesis were compiled using this batch file. 

File ACOMP.BAT 

PATH=\ms_c\wrk;\ms_c\tmp;\ms_c\src;\ms_c\bin;\ms_c\lib;\ms_c\include;c:\dos 
SET INCLUDE =\MS_C\INCLUDE 
SET LIB=\ms_c\lib 
SET TMP=\ms_c\tmp 

cl /c /Fs /Ze /GO /Alnd %1 .c 

IF NOT ERRORLEVEL 1 link /NOD /E %1,„rTEXPC.UB support.lib mlibce.lib; 
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